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ABSTRACT 


Multiple  linear-regression  prediction  equations  have  been  derived  for  the  winter¬ 
time  12- ,  24-,  and  36-hr  displacements  and  changes  in  central  pressure  of  cyclones  and 
anticyclones  over  Europe  and  cyclones  over  Asia.  An  inq>ortant  feature  of  the  statistical 
prediction  technique  used  in  developing  the  equations  is  that  the  coordinates  of  measure¬ 
ment  accompany  the  cyclone  or  anticyclone.  The  equations  were  tested  on  sizable  inde¬ 
pendent  data  sanq>les  and  found  to  be  stable;  limited  coixq>arative-accuracy  tests  indicate 
that  their  results  are  conq>etitive  with  those  obtained  by  e;q>erienced  subjective  forecasters. 
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1.0  INTRODUCTION 


Predicting  the  behavior  of  significant  features  in  the  distribution  of  surface  pres¬ 
sure,  such  as  cyclmes,  anticyclones,  fronts,  ridges  (high  pressure),  and  troughs  (low 
pressure),  is  an  important  step  in  the  preparation  of  a  weather  forecast.  Gross  sig¬ 
nificant  features  like  these  have  drawn  attention  because  their  initial  configuration  and 
their  subsequent  movement  and  intensification  (or  dissipation)  permit  a  reasonably 
accurate  sequential  specificaticm  of  the  weather  at  stations  along  their  paths.  These 
relationships  have  long  been  regarded  as  primarily  qualitative.  Their  quantitative  study, 
which  is  urgently  needed,  can  lead  to  new  characteristics  that  should  be  predicted  or  to 
new  significant  features  that  are  rich  in  weather-predictive  information. 

Formulating  the  problem  of  predicting  the  distribution  of  surface  pressure  in  terms 
of  the  behavior  of  significant  features  permits  two  important  departures  from  the  custom¬ 
ary  statistical  formulation,  in  which  predictors  for  a  given  gric^int  pressure  or  array 
of  gridpoint  pressures  are  defined  as  observations  in  a  fixed-coordinate  system  [6] . 

First,  it  permits  the  minimization  of  the  errors  (in  a  least-squares  sense)  of  more 
meaningful  parameters  than  simply  the  pressure  at  a  gridpoint.  Second,  it  permits  the 
introduction  of  a  type  of  nonlinearity  in  the  predictor-predictand  relationships.  Consider, 
for  example,  the  extension  of  this  iqiproach  to  the  derivation  of  equations  for  forecasting 
the  behavior  of  a  large  number  of  significant  features.  The  forecasts  can  be  combined 
by  methods  analogous  to  those  used  in  objective  weather  analysis  to  yield  a  prediction  of 
the  surface  pressure  at  an  array  of  gric^oints.  The  predicted  pressure  at  any  gridpoint 
is  a  function  of  the  synoptic  situation;  that  is,  predictors  affecting  the  pressure  at  the 
gridpoint  stem  from  the  location  and  type  of  significant  features  present  at  the  initial 
time  of  the  forecast  period. 

In  the  work  described  here,  the  cyclone  (defined  as  having  at  least  one  closed  isobar) 
was  chosen  as  the  primary  feature  for  study.  Partly  because  of  its  association  with 
cloudiness  and  precipitation,  the  cyclone  has  been  under  investigation  by  meteorologists 
since  the  advent  of  synoptic  maps.  The  anticyclone  (likewise  defined  as  having  at  least 
one  closed  isobar)  was  also  studied.  Earlier  work  on  the  problem  of  predicting  the  move¬ 
ment  and  intensification  of  cyclones  over  the  eastern  United  States  by  statistical  methods 
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{ 12]  proved  promising  enough  to  warrant  the  extension  of  such  an  i4>proach  to  other 
geographical  areas.  Accordingly.  Europe  and  eastern  Asia  were  selected  as  regions  for 
technique  development. 
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2.0  DESIGN  OF  THE  STUDIES 


2.1  Areas  Studied 

Regression  equations  for  predicting  the  wintertinne  (November-March)  behavior 
of  cyclones  were  derived  separately  for  Europe  and  eastern  Asia  and  for  each  of  the 
zones  defined  in  Fig.  2-1.  Anticyclones  were  examined  also,  and  equations  were  derived 
for  Europe;  the  "winter  reason"  was  extended  to  include  October  and  April  so  that  there 
would  be  enough  cases  to  test. 


Fig.  2-1.  Areas  for  which  cyclones  and  anticyclones  were  predicted.  Area  I  (Europe) 
is  divided  into  two  zones  (northern  and  southern)  by  latitude  45“N.  Area  II  (Asia)  is  divided 
into  four  zones  (northwestern,  northeastern,  southwestern,  and  southeastern)  by  latitude 
40“ N  and  longitude  14 5° E. 
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TABLE  2-1 

SELGCTEO  MINTER  CYCLONES  AND  ANTICYCLONES  OVER  EUROPE  AND  ASIA, 
1955-1960 


ArM* 

Zona 

No.  of  cyelonai 

No.  of  anti eye lonas 

1955-1999t 

1959-1960* 

1955-1960 

1955-19591 

1959-1960* 

1955-1960 

N 

402 

121 

523 

- 

- 

- 

Europ«|l 

S 

518 

82 

4x 

- 

- 

- 

Both 

720 

203 

•  925 

656 

106 

764 

NH 

266 

72 

538 

- 

- 

- 

NE 

435 

lie 

551 

- 

- 

- 

Aslal 

SN 

196 

59 

255 

- 

- 

- 

SE 

196 

57 

255 

- 

- 

- 

Al  1 

1091 

266 

1557 

- 

- 

•Sm  F^g.  2-1. 

tOapandanf  sanpU. 

tindapandant  tamp  la. 

loividing  llna  batoaan  zonas  It  1l5*N. 

lOlvIdIng  I Inai  batwaan  zonas  are  4o*N  and  l4;*E. 


2.2  Selectioi.  of  Cases 

Table  2- 1  shows  the  number  of  cyclones  and  anticyclones  selected  for  each  of  the 
zones  of  Fig.  2-1.  The  cyclones  were  selected  by  examining  all  the  0000-  and  1200-GCT 
surface  charts**  for  the  winters  (November— March)  of  1955-1956  through  1959-1960. 

A  cyclone  was  accepted  if  it  retained  its  identity  for  at  least  36  hours.  European  anti¬ 
cyclones  were  selected  in  the  same  way  except  that  the  "winter  season"  was  extended  to 
include  October  and  April  to  provide  a  large  enough  sample. 

The  samples  of  cyclones  and  anticyclones  selected  for  the  winters  of  1955-1956 
through  1958-1959  were  designated  as  the  dependent  samples  and  were  used  in  deriving 
the  equations.  The  samples  for  the  winter  of  1959-1960  were  designated  as  the  inde¬ 
pendent  samples  and  were  used  to  test  the  equations. 


**Microfilm  copies  of  surface  maps,  analyzed  by  the  National  Weather  Analysis 
Center,  were  used.  From  these  it  was  possible  to  construct  maps  showing  individual 
cyclone  and  anticyclone  tracks  at  12-hr  intervals,  which  served  as  the  basis  for  selec¬ 
tion. 
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TABLE  2-2 
PREDICTANOS 


Class 

Symbol 

Unit  of 
measurement 

Northward  displacement 

A 

N 

deg.  lat. 

Eastward  displacement 

E 

deg.  lat. 

Change  in  central  pressure 

5 

mb 

2.3  Predictands 

Latitudinal  displacement,  longitudinal  displacement,  and  change  In  central  pres¬ 
sure,  for  12,  24,  and  36  hr,  were  chosen  as  predictands  (Table  2-2)  for  both  the  cyclone 
and  anticyclone  equations.  The  same  National  Weather  Analysis  Center  mtqis  used  to 
select  the  cases  (see  Sec.  2.2)  were  used  as  the  source  of  predictand  data. 

2.4  Predictors  Considered 

System  433L  hemispheric-data  tapes  [8]  were  used  as  the  basic  source  of  predictor 
data.  Special  preprocessing  programs  (see  App.  A)  automatically  derived  grid^oint  arrays 
of  the  various  pressure,  height,  and  thickness  data  for  each  cyclone  and  anticyclone  In  the 
developmental  sample,  yielding  altogether  six  grid  arrays  of  221  points  each,  or  about 
1300  potential  predictors. 

To  comply  with  the  maximum  of  200  possible  predictors  acceptable  by  the  screening- 
regression  program  (see  Sec.  2.6),  a  preliminary  reduction  in  the  number  of  predictors 
was  made  subjectively  by  increasing  the  spacing  between  gri(4>oint8.  Due  to  the  spatial 
redundancy  of  the  meteorolc^ical  information  under  consideration,  this  hypothetically  lost 
very  little  of  the  predictive  information  contained  in  the  original  set  of  data. 

Table  2-3  lists  the  possible  predictors  that  remained.  There  are  two  categories: 
one  includes  both  surface  and  iqjper-air  predictors  in  combination,  and  the  other  includes 
surface  predictors  only.  No  attempt  was  made  to  incorporate  so-called  derived  predictors 
(vorticity,  vorticity  advectlon,  etc.)  because  a  previous  feasibility  test  involving  such  pre¬ 
dictors  1 12)  did  not  demonstrate  significant  improvement  over  simply  using  point-value 


TABLE  2-3 

POSSIBLE  PREDICTORS 


Class 

Symbol 

Unit  of 
measurement 

No.  of  predictors 

Including 

upper-air 

data 

Excluding 

upper-air 

data 

Ssa-laval  prassura 

P 

mb 

55 

84 

12-hr  saa-laval  prassura  change 

AP 

mb 

J4 

84 

300-nib  height 

Z 

10  ft 

4i 

- 

12-hr  300-mb  height  change 

10  ft 

4o 

- 

1000-300-mb  thickness 

H 

10  ft 

24 

- 

12-hr  1000-300-mb  thickness  change 

10  ft 

24 

- 

Present  latitude  ot  cyclone 

deg.  lat. 

1 

1 

Present  longitude  of  cyclone 

\ 

_ _ _ _ i 

deg.  long.* 

1 

1 

Total 

200 

170 

•Positive  tor  W  long.,  negative  tor  E  long. 


predictors,  although  only  a  thorough  investigation  could  support  any  contention  that  derived 
predictors  contain  no  independent  predictive  information. 

2.5  The  Moving-coordinate  Grid 

The  grid  for  extracting  predictor  infor.  ation  accompanies  the  feature  as  it  moves, 
so  variables  are  measured  at  constant  positions  relative  to  its  center.  The  grid  is  shown 
in  Fig.  2-2.  The  gridpoint  defined  by  the  (k,f)- location  (10,  5)  is  placed  at  the  center  of 
the  feature,  and  the  grid  is  oriented  so  that  the  line  k  =  10  coincides  with  the  meridian 
passing  through  the  center  of  the  feature.  In  practice,  grid  placement  and  data  tabulation 
are  done  by  computer  programs,  and  "analyzed  maps"  are  on  magnetic  ts^e.  On  a  polar 
stereogrs^hic  projection  with  standard  parallel  at  60’N,  the  17  x  13  array  forms  a  set  of 
evenly  spaced  points.  On  a  map  scale  of  1:15,000,000,  the  grid  interval  is  1  in.  This  is 
the  same  interval  as  in  the  JNWP  grid  [2,  Fig.  1).  At  60'N,  where  the  scale  is  true,  one 
grid  interval  equals  381  km.  The  221  points  defined  by  this  grid  system  were  the  ones 
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Fig.  2-2.  Moving-coordinate  grid  overlay.  Gridpoint  (10,  5) 

in  the  (k,  i)-array  is  always  positioned  at  the  center  of  the  cyclone 
or  anticyclone,  and  the  grid  is  always  oriented  so  that  the  line 
k  =  10  coincides  with  the  meridian  passing  through  the  center  of 
the  feature.  One  grid  interval  equals  one  JNWP  grid  interval 
(381  km  at  60“  N). 

used  for  predictor  tabulation  by  superimposing  the  grid  on  analyzed  maps  of  sea-level 
pressure  and  500-mb  height. 

An  important  preliminary  step  in  the  preparation  of  the  data  for  the  prediction 
experiments  was  a  systematic  machine-man  search  for  and  correction  of  gross  errors. 
Suspicious  gridpoint  values  of  parameters  were  identified  by  a  series  of  tests  performed 
by  the  IBM  7090,  and  listings  of  these  values  were  compared  with  microfilm  copies  of 
manuscript  maps;  where  necessary,  corrections  were  made  to  the  data  filed  on  magnetic 
tapes.  The  details  of  the  computer  programs,  data-handling  procedures,  and  error¬ 
checking  procedures  are  given  in  .Vppendix  A. 

2.G  Screening-regression  Technique 

The  screening  procedure  suggested  by  Bryan  [1]  and  developed  for  the  IBM  704 
electronic  computer  by  Miller  |5,  7)  was  used  to  screen  the  possible  predictors  identified 
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in  Table  2-3.  One  who  designa  a  statistical  prediction  eiqperiment  likes  invariably  to 
consider  predictors  deemed  important  on  the  basis  of  previous  theoretical,  synoptic, 
and  empirical  work  (in  this  case,  all  those  listed  in  Tabie  2-3),  but,  as  Lorenz  [4]  points 
out,  a  prediction  equation  should  contain  few  predictors  in  comparison  with  the  size  of 
the  developmental  sanq)la;  if  there  are  too  many,  a  relationship  that  fits  the  sample  used 
to  establish  it  is  likely  to  fail  when  tqiplied  to  a  new  sample.  The  object  of  the  screening 
procedure  is  to  select  from  a  set  of  possible  predictors  the  subset  that  contributes  most 
significantly  and  independently  to  reducing  the  variance  of  the  predictand. 

From  an  array  of  possible  predictors,  the  screening  procedure  selects  first  the 
one  that  has  the  highest  linear  correlation  with  the  predictand  in  question.  This  predictor 
is  then  held  constant,  and  partial-correlation  coefficients  between  the  predictand  and  each 
of  the  remaining  predictors  are  examined;  the  predictor  now  associated  with  the  highest 
coefficient  is  the  second  one  selected.  Additional  predictors  are  chosen  similarly  until 
a  selected  predictor  fails  to  explain  a  significant  additional  percentage  of  the  remaining 
variance  of  the  predictand. 

The  criterion  of  significance  as  applied  to  the  screening  procedure  is  not  clear  cut 
since  the  usual  F-test  methods  (e.g.,  [10])  are  not  ^plicable  [7].  If  a  predictor  is  chosen 
at  random  from  a  group  of  predictors,  an  F-test  is  usually  taken  at  the  95%  level;  this 
allows  for  a  l-in-20  chance  of  considering  the  predictor  significant  when  in  fact  it  is  not. 
Because  the  screening  procedure  does  not  select  predictors  randomly,  a  more  severe  test 
is  needed  to  specify  a  l-in-20  chance.  For  his  screening  procedure.  Miller  suggested 
[7]  that  the  critical  F-value  be  a  function  of  the  number  of  possible  predictors.  The  F-test 
was  used  in  this  form  in  these  e;q>eriments. 


8 


3.0  RESULTS 


Preparatory  to  deriving  prediction  equation!  for  cyclonea  or  smticyclonea.  the 
climatological  characteristics  of  tiiese  features  and  the  atmospheric  variables  believed 
to  be  related  to  them  were  conq>uted  as  an  aid  for  suggesting  possible  predictors  to  be 
used  in  the  screening-regression  program.  Also  of  interest  was  a  determination  of  the 
difference  in  characteristics  for  different  regions  of  the  same  area  to  assess  the  desira¬ 
bility  of  stratifying  the  data  sample.  Accordingly,  the  means  and  standard  deviations  of 
all  predictands  and  possible  predictors  were  computed  both  for  entire  areu  (unstratified) 
and  for  regions  within  areas  (stratified).  Some  of  these  characteristics  are  discussed  in 
the  sections  that  follow,  in  addition  to  the  results  obtained  by  the  regression  analysis. 

The  prediction  equations  are  given  in  Appendix  B. 

3. 1  European  Cyclones 

Europe  (area  I  of  Fig.  2-1)  was  divided  by  latitude  45*N  into  two  aones  on  the  basis 
of  charts  of  principal  tracks  of  cyclones  and  their  frequencies  described  by  Klein  [3], 
which  indicate  separate  maxima  of  cyclone  activity  and  different  behavior  over  northern 
Europe  and  in  the  Mediterranean  area.  The  noain  purpose  of  the  stratification  was  to 
attempt  to  account  for  differences  in  predictor-predictand  relationships  that  may  be  present 
between  one  geogr^)hical  zone  and  another.  The  need  for  stratification  is  not  always 
clear  cut.  For  example,  it  is  possible  that  one  or  more  of  the  predictors  themselves  may 
be  able  to  eiqpress  such  differences,  making  stratification  unnecessary. 

3.1.1  Synoptic  Climatology 

Table  3- 1  contains  the  means  and  standard  deviations  of  the  northward  and  eastward 
displacements  amd  changes  in  central  pressure  of  720  cyclones  for  12,  24  and  36  hr.  This 
information  is  also  shown  grsqjhically  in  Fig.  3-1.  The  mean  tracks  were  constructed  by 
simply  connecting  the  mean  displacements  for  the  three  time  periods  of  interest.  Note 
that  the  mean  track  for  both  northern-  and  southern-zone  cyclones  is  toward  the  east- 
northeast  but  that  the  speed  of  the  northern  cyclones  (about  15  knots)  is  slightly  greater 
than  that  of  the  southern  cyclones  (about  12  knots).  Both  sanqrles  exhibit  filling  character¬ 
istics  in  contrast  to  the  deepening  tendencies  of  cyclones  in,  say,  the  eastern  United  States. 
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TAas.E  3-1 

CmRACTERISTICS  OF  WINTER  CYCLONES  OVER  EUROPE. 
1935-1959  (OEPENOENT  SAMPLE) 


Zone 

Forecast 

Interval, 

hr 

Observed 
nort heard 
displacement, 
deg.  lat. 

Observed 
essteard 
displacement, 
deg.  1st. 

Observed 
change  In 
central  pressure, 
mb 

Mean 

Std.  dev. 

Maan 

Std.  dev. 

Mean* 

Std.  dev. 

12 

0.57 

2.39 

.3.22 

2.49 

-0.09 

5.82 

N 

24 

i.oe 

4.00 

.6.02 

4.13 

1.51 

8.70 

36 

1.96 

5.27 

-8.71 

5.68 

3.74 

10.08 

12 

0.4o 

2.09 

-2.38 

2.13 

-0.24 

3.99 

S 

24 

0,85 

3.15 

-4.77 

5.48 

0.66 

5.46 

56 

1.32 

4.o4 

-6.94 

4.76 

2.14 

6.62 

12 

0.38 

2.26 

-2.85 

2.58 

-0.16 

5.10 

Both 

24 

0.98 

5.65 

-5.47 

5.91 

1.14 

7.46 

56 

1.69 

4.78 

-7.95 

5.56 

3.03 

9.21 

•fl«gative  values  represent  deepening. 


a 

CD 

T, 

N 

i 

JB 

t: 

o 

Z 


Fig.  3-1.  Mean  tracks  of  winter  cyclones  in  Europe  by  zone,  1955- 
1959  (dependent  sample).  oi2-hr  displacement;  A24-hr  displacement; 
•36-hr  displicement.  Value  adjacent  to  symbol  refers  to  mean  change 
in  central  pressure  (millibars). 
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S  (0:39.1  9^:-l3.8) 

Pig.  3-2.  Mean  maps  of  sea-level  pressure  for  winter  cyclones  in 
Europe,  1936-1959  (dependent  sample).  Isobars  are  labeled  in  millibars. 


5  (0:39.1  •A:-I3.8) 


Fig.  3-3.  Mean  maps  of  12-hr  pressure  change  for  winter  cyclones  in 
Europe,  1955-1939  (dependent  sample).  Issdlobars  are  labeled  in  millibarb. 
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5  (0:39.1  7^:-l3.8) 

Fig.  3-4.  Mean  maps  of  500-mb  height  for  winter  cyclones  in  Europe, 
1955-1959  (dependent  sample).  Isohypses  are  labeled  in  tens  of  feet. 
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Fig.  3-5.  Mean  maps  of  12-hr  500-mb  height  change  for  winter  cyclones 
in  Europe,  1955-1959  (dependent  sample).  Isallohypses  are  labeled  in  tens 
of  feet. 
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Fig.  3-6.  Mean  maps  of  1000-500-mb  thickness  for  winter  cyclones  in 
Europe,  1955-1959  (dependent  sample).  Isopachs  are  labeled  in  tens  of  feet. 
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Fig.  3-7.  Mean  maps  of  12-hr  1000-500-mb  thickness  change  for  winter 
cyclones  in  Europe.  1955-1959  (dependent  sample).  Isallc^achs  are  labeled 
in  tens  of  feet. 
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Although  the  cyclone  tracke  in  the  two  zones  seem  similar,  there  is  a  sizable  difference 
in  variability  about  the  means,  as  indicated  by  the  standard  deviations  in  Table  3-1.  For 
every  time  period,  the  standard  deviations  for  the  northern  zone  are  greater  than  for  the 
southern.  Other  interesting  characteristics  are  shown  in  Figs.  3-2  through  3-7,  in  which 
the  means  of  all  the  possible  predictors  have  been  plotted  and  analyzed.  The  resulting 
charts  depict  the  mean  distributions  of  pressure,  height,  and  thickness  and  their  12-hr 
changes  relative  to  the  cyclone  center.  Comparison  of  Figs.  3-2(a)  and  (b)  shows  that 
the  mean  northern-zone  cyclone  (984  mb)  is  about  17  mb  deeper  than  the  mean  southern- 
zone  cyclone  (1001  mb)  and  is  much  larger  in  extent.  This  is  characteristic  of  most  of 
the  other  features;  i.e.,  the  northern  zone  shows  pressure-  and  height-change  centers  of 
greater  magnitude,  and  higher-amplitude  500-mb  troughs. 

3.1.2  Results  of  Using  Both  Surface  and  Upper-air  Predictors 

Regression  equations  were  derived  for  the  two  zones  separately  and  for  the  entire 
area.  The  possible  predictors  that  were  screened  are  listed  in  Table  2-3.  The  predict- 
ands  are  northward  displacement,  eastward  displacement,  and  change  in  central  pressure, 
for  12,  24,  and  36  hr. 

The  results  of  the  screening  regression  on  dependent  data  are  presented  in  Tables 
3-2  and  3-3.  Table  3-2  lists  the  predictors  in  the  order  of  their  selection  by  the  screening 
procedure  and  the  percentage  of  the  total  variance  of  the  predictand  e^lained  by  each  for 
northern-zone  cyclones,  southern-zone  cyclones,  and  the  complete  unstratified  sample. 
Table  3-3  summarizes  the  initial  and  residual  standard  deviations,  along  with  the  percent 
reductions.  The  numbers  in  parentheses  accompanying  the  predictors  in  Table  3-2  refer 
to  the  (k,  i)-predictor  locations  in  the  grid  system  shown  in  Fig.  2-2. 

From  Table  3-2(a),  it  can  be  seen  that  the  six  displacement  predictands  for  the 
nothern  zone  (two  fcr  each  time  period)  counted  heavily  on  iq)per-air  information.  For 
the  northward  displacement,  Z(15,  5)  was  the  first  predictor  selected  for  all  three  time 
periods.  Z(15,  5)  is  the  500-mb  height  about  six  grid  intervals  east  of  the  cyclone.  It  is 
correlated  with  the  northward  component  in  such  a  way  as  to  favor  large  northward  dis¬ 
placements  when  the  heights  there  are  high.  This  seems  consistent  with  synoptic  reasoning; 
i.e.,  a  large-amplitude  ridge  east  of  a  cyclone  is  conducive  to  a  sizable  northward  displace- 
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TABIE  5-2 

PREDICTORS  SELECTED  BV  SCREENING  REGRESSION 
FOR  M INTER  CVCLONES  OVER  EUROPE* 

(a)  Northern  zona 


Foraeatt 

Ord«r  of 

N 

E 

0 

Intarval t 

% 

hr 

sal  act  Ion 

% 

% 

Predictor 

rad  a 

Pradlctor 

rad. 

Predictor 

rad. 

let 

T(15,5) 

25.3 

2(  9,1) 

26.8 

fi.P(  9,5) 

26.6 

?nd 

2( 

15.0 

2(15,9) 

10,2 

P(10,5) 

9.1 

5rd 

P(  5,1) 

5.8 

'SP(  9,5) 

5.1 

PCi1,5) 

2.9 

Itth 

2(13,5) 

4.4 

2(15,1) 

2.1 

Pf15,l) 

5.1 

12 

5th 

- 

- 

2(11,7) 

2.0 

6P(  7,1) 

2.5 

6th 

- 

- 

2(  9,5) 

2.5 

62(  9,5) 

1.9 

- 

- 

- 

- 

A 

2.1 

8th 

- 

- 

- 

■ 

2(  9.7) 

2.5 

Total 

48.5 

48.7 

50.5 

1st 

2(15,5) 

25.5 

2(  9,1) 

26.7 

07(  9,5) 

20.0 

2nd 

2(  9,7) 

15.6 

7(15,9) 

17.1 

P(10,5) 

15.2 

Srd 

7(  5,1) 

5.5 

2(15,1) 

4a1 

OH(  9,5) 

6.4 

2l| 

l*th 

7(11,1) 

4.9 

^2(  9,5) 

2.7 

A 

2.7 

•> 

2.6 

2(15,7) 

1.6 

2(  9,7) 

5.6 

6th 

P(15,7) 

1.5 

P(I0,5) 

1.2 

- 

- 

7th 

H(  9,1 ) 

1.5 

2(5,11) 

1.7 

- 

Tola  1 

56.7 

55.1 

45.9 

1st 

2(15,5) 

20.4 

2(  9.1) 

24.4 

7(  9,3) 

19.3 

2rtd 

7(  9,7) 

>7.1 

2(15,9) 

OP(  9,5) 

11.6 

5r<l 

'2(15,5) 

5.8 

2(15,1) 

P(io,5) 

9.6 

'56 

l|th 

7(  7.1) 

2.9 

'>2(  9.5) 

m 

H(  9,7) 

2.3 

5th 

.  t 

2.8 

7(15,7) 

1.7 

2(15,1) 

2,8 

Stt. 

2(15,1) 

2.4 

- 

P(li,5) 

2.7 

Total 

||9.4 

50..? 

48.5 

Inc.'jd  ng  u^per-oir  drlt 


n 

ii 


(b)  teutharn  ten* 


P^rwwt 

Intarval, 

hr 

fl 

E 

_ 6 _ 1 

Hi 

H 

fradletor 

H 

fradictor 

H| 

lat 

7,7) 

7(11,9) 

gg 

6P(  9,5) 

HDji 

2a4 

T(I3,5) 

M 

^(  5,7) 

^(10,5) 

13' 

Vb 

7{  9,7) 

■9 

7(  7,1) 

■n 

P(11,7) 

■3 

Total 

*k.7 

21.1 

26.k  1 

lat 

K  7,5) 

11.7 

7(11,9) 

2(10,5) 

16.8 

2n4 

T(I9,9) 

6.2 

^(  5,7) 

n 

02(  9,5) 

15.6 

2k 

Jrb 

n  9,9) 

12.2 

7(  9,5) 

k.i 

2(11,7) 

kth 

AT(  7,9) 

■g 

- 

- 

» 

M 

5th 

♦ 

M 

- 

- 

H(  9,9) 

M 

Total 

55.1 

•77.5 

kk.3 

ut 

T(  9,9) 

7(15,9) 

16.5 

2(10,5) 

2nd 

7(15,5) 

P(  5,7) 

10.7 

7(11,9) 

Vd 

H  5,5) 

m 

• 

- 

\ 

kth 

- 

m 

- 

• 

02(  9,5) 

5th 

• 

• 

- 

0P(  5,5) 

1  Total 

55.9 

77.7 

)l6.6  I 

i 

) 


I 
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(c)  Both  zona* 


Forarait 

Intarvati 

hr 

Ordar  o« 
aataetion 

N 

1 

6 

Pradictor 

im 

Pradictor 

% 

Pradictor 

% 

rad. 

Ut 

P(  7,1) 

P(15,9) 

m 

6P(  9,5) 

mmm 

2fld 

P(15,3) 

P(11,1) 

ig 

P(10,5) 

■1 

Jrd 

P(  9.7) 

s.li 

P(  5,7) 

7(  9,7) 

■i 

kfh 

7(15,5) 

9.0 

&P(  7.5) 

5.4 

X 

5.4 

5th 

M(  7,7) 

li.k 

7(  9,1) 

2.1 

P(11,5) 

1.8 

12 

6th 

P(ll,3) 

1.8 

7(11,7) 

5.2 

p;i9,3) 

2.0 

7th 

- 

- 

P(10,3) 

2.6 

07(  9,5) 

1.4 

8th 

- 

• 

P(’l,5) 

1.9 

- 

- 

9th 

- 

- 

H(13,3) 

1.5 

- 

- 

tOth 

- 

- 

7(  5,9) 

1.5 

- 

- 

11th 

II 

- 

- 

P(  9,5) 

1.2 

. 

Total 

39.7 

43.4 

4i.o 

lat 

P(  7,1) 

10.9 

P(15,9) 

16.8 

.  aP(  9,5) 

15.6 

2nd 

P(15,3) 

7.5 

.  P(11,1) 

7.7 

P(10,5) 

12.5 

yd 

P(  9,9) 

5.0 

P(  5,7) 

7(  9,7) 

6.0 

llth 

7(15,5) 

7.6 

7(  9,1) 

\ 

5.6 

2>l 

yth 

7(  9,7) 

8.2 

7(13,9) 

^9 

0P(11,5) 

2.1 

6th 

7(11.3) 

3.3 

07(  7,5) 

1.8 

P(11,5) 

7th 

A7(  9,9) 

1.7 

7(11,9) 

0a6 

P(13,1) 

0th 

M(  7,3) 

1.6 

•  h(15,3) 

laO 

a7{  9,5) 

9th 

P(ll,7) 

7(15,7) 

1.7 

OP(ll,9) 

10th 

P(15,7) 

n 

P(10,5) 

2.1 

• 

Total 

<•7.9 

45.9 

46.4 

lit 

P(  7,1) 

9.3 

P(15,9) 

18.9 

P(10,5) 

16.1 

2nd 

P(  9,9) 

7.3 

P(il,l) 

7.2 

iP(  9,5) 

i4.5 

3rd 

7(15,5) 

7.7 

P(  5,7) 

4.9 

7(15,9) 

6.1 

llth 

7(  9,7) 

10.6 

7(  9,1) 

1.9 

2.7 

ye 

5th 

P(15,l) 

5.5 

7(13,9) 

5.0 

7(  9,7) 

3.1 

6th 

a7(  9,9) 

1.5 

OP(  9,5) 

1.5 

A7(ll,3) 

1.3 

7th 

^17(13,3) 

l.li 

7(15,5) 

laO 

H(11,3) 

1.5 

0th 

- 

• 

P(  9,3) 

0.7 

07(11,7) 

1.0 

9th 

- 

- 

P(15,5) 

1.1 

- 

- 

Totol 


<•6.3 


meat.  Z(9. 1),  the  500-xnb  height  about  four  grid  intervals  south-southwest  of  the  cyclone, 
was  selected  first  for  all  three  time  periods  for  the  eastward  component.  Its  correlation 
with  the  eastward  component  is  such  that  low  heights  in  this  region  are  associated  with 
small  eastward  displacements.  Low  heights  in  that  location  imply  high-an4>litude  sys¬ 
tems  at  500  mb  or  lack  of  a  str(»g  sonal  flow.  The  central  pressure's  lead  predictors  of 
P(9,  S)  for  12  hr  and  Z(9. 5)  for  24  hr  correlate  12-hr  pressure  and  height  falls  to  the 
immediate  west  of  the  cyclone  with  deepening. 

In  the  southern  none  (Table  3-2(b)] ,  there  is  a  greater  reliance  placed  on  surface 
predictors.  This  is  carried  over  even  to  the  unstratified  equaticms  [Table  3-2(c)],  where 
the  first  (and  even  the  sectmd)  predictor  selected  is  a  surface  variable  in  all  nine  instances. 
For  the  northward  component,  P(7,l),  the  sea-level  pressure  five  grid  intervals  southwest 
of  the  cyclone,  correlates  high  pressure  there  with  large  northward  displacements. 

P(15, 9),  the  lead  predictor  for  eastward  displacement,  is  located  about  six  grid  intervals 
northeast  of  the  cyclone,  and  the  relationship  is  such  that  high  pressure  there  retards 
eastward  cyclone  displacement.  For  central  pressure,  P(9,  5)  is  the  same  predictor  as  in 
the  northern  zone.  P(10, 5)  is  the  initial  central  pressure  itself  and  is  correlated  in  such 
a  way  that  little  additional  deepening  takes  place  when  the  cyclone  is  initially  de^. 

Of  particular  interest  is  the  conq>arison  between  the  stratified  and  unstratified 
regression  analyses.  Such  a  con^arison  is  shown  in  Table  3-4.  The  vector  column  refers 
to  the  rms  vector  error,  which  is  the  square  root  of  the  sum  of  the  squares  of  the  two 
component  rms  errors.  Note  that  the  vector  errors  for  the  unstratified  equations  are 
smaller  for  all  three  time  periods  for  both  the  northern  .time  and  the  southern  zone.  For 
central  pressure,  the  unstratified  equations  yield  smaller  rms  errors  exc^t  for  the  36-hr 
southern- zone  predictions. 

The  failure  of  the  stratified  equations  to  inqirove  the  results  would  seem  to  indicate 
that,  at  least  for  this  problem,  there  is  no  benefit  in  processing  the  cases  separately 
according  to  the  latitude  of  the  cyclone.  One  factor  to  consider  is  that  720  cases  were 
\ised  to  develop  the  unstratified  equations,  whereas  only  about  half  that  number  were  used 
to  develop  the  stratified  equations  for  each  zone.  The  smaller  number  of  cases  could 
result  in  less  stable  coefficients.  If  this  is  true,  then  much  care  needs  to  be  taken  in 
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TABLE  3-5 

CUMULATIVE  FREQUENCY  OF  ERROR  IN  24.hr 
PREDICTIONS  OF  CYCLONES  OVER  EUROPE, 
1959-1960* 


(a)  Total  dlaplaeamant 


Error  c , 
dag.  iat. 

Parcantaga  of  casas 

Ragrassiont 
(2^  casas) 

Dunstable 
(784  eases) 

Offenbach 
(1195  cases) 

Moscow 
(1126  cases) 

mam 

8 

6 

6 

8 

23 

19 

18 

23 

43 

37 

35 

43 

61 

53 

51 

61 

0  ^  e  s5 

75 

64 

61 

72 

0  s  e  26 

87 

74 

72 

81 

OSes? 

90 

85 

84 

90 

(b)  Change  in  central 

pressure 

Error  c , 
mb 

Percentage  of  cases 

Ragrassiont  . 
(2(9  casas) 

Dunstable 
(784  cases) 

Offenbach 
(1195  cases) 

Moscow 
(1126  cases) 

25 

22 

23 

28 

mm 

61 

54 

58 

64 

^B9 

79 

75 

77 

83 

0  S  €  SIO.5 

92 

84 

87 

92 

0  s «  ^13.5 

98 

91 

92 

96 

0  s  e  SI6.3 

99 

96 

96 

98 

0  s  e  SI9.3 

99 

96 

98 

99 

*Ragraasion  pradietions  art  for  Nov.  1959  through  Mar.  I960.  Ounstabla, 
Offanbach,  and  Moscow  pradietions  ara  for  Jun.  1959  through  Mar.  i960. 
tUnstratIf lad. 


24 


determining  oriterU  for  etretlfloation.  pertiouUrly  if  the  inveetlgetor  has  a  limited 
amoimt  oi  data  at  hia  diapoaal. 

Table  3-4  alao  oootaina  a  ooxq>ariaon  of  the  regreaaion  prediotiona  with  olimatology. 
The  olimatology  uaed  waa  baaed  on  the  computed  meana  (TaUe  S-1)  of  the  prediotande  for 
the  dependent  Bang>le.  Theae  were  uaed  aa  prediotiona  for  the  203  oaaea  of  the  independent 
aan^le.  The  rma  errora  for  olimatology  are  aubatantially  larger  than  for  regreaaion  for 
all  time  perioda. 

To  determine  the  uaefUlneea  of  a  prediction  technique,  it  ia  deairable  to  teat  ita 
aklll  againat  aome  other  teohniquo^referably  the  current  operational  ayatem  for  which 
the  new  technique  ie  hoped  to  aerve  aa  an  aid.  One  auch  compariaon  would  be  with  prog- 
noaea  iaaued  by  a  foreoaat  center.  Thia  type  of  teat  la  beii^  performed  by  a  U8AF  Air 
Weather  Service  (AW8)  unit  in  Eur(^,  but  reaulta  are  not  currently  available.  Some  indi¬ 
cation  of  the  performance  of  the  regreaaion  techniquea  can  be  obtained,  though,  from  a 
conqpariaon  of  aome  verlficationa  of  prognoatio  charta  prepared  by  Soherhag  [9] .  In  hia 
report.  Soherhag  preaenta  error  diatributiona  of  24-hr  cyclone-poaition  and  central- 
preaaure  foreoaata  for  three  forecast  centera:  Dunatable.  Offenbach,  and  Moscow.  These 
results  are  shown  in  Table  3-5.  with  the  203  independent  regression  predictions  (unstrati¬ 
fied)  for  comparison.  Although  the  data  samples  are  not  exactly  the  same.  Table  3-5 
indicates  that  the  regression  technique  at  least  maintained  the  same  level  of  skill  as  the 
prognoses  of  the  Dunstable.  Offenbach,  and  Moscow  forecast  centers. 

3.1.3  Results  of  Using  Surface  Predictors  Only 

There  are  limitations  to  the  application  of  this  prediction  technique.  For  exanq>le. 
since  SOO-mb-height  data  are  received  only  every  12  hr  and  the  method  requires  iiq>ut  of 
such  information,  it  is  not  possible  to  apply  the  regression  equations,  (^eratlonally,  at  0600 
and  1800  GOT.  Also,  in  many  instances,  there  is  a  time  lag  of  at  least  an  hour  between  the 
receipts  of  surface  and  unper-air  information  for  the  aiqpUcable  times  of  0000  and  1200  OCT. 
One  possible  way  to  overcome  these  problems  is  to  use  a  prognosis  of  the  500-mb  data. 
Another  approach  is  to  derive  r^ession  equations  that  use  only  surface  data  as  predictors. 
A  regression  analysis  that  excludes  iq>per-air  data,  if  successful,  would  provido  a  useful 
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table  ]-6 

MEOICTOM  SELECTED  BT  SCREENIN8  MBSIIESSION 
FOB  WINTER  CVCIONES  (WEN  EUROPE* 


POTMMt 

Ordor  of 
•aiactlon 

N 

E 

6 

1 V9 1 1 
hr 

Prodlctor 

lEP 

Prodletor 

% 

r»dp 

Prtdlctor 

% 

rod. 

Ut 

P(  7,1) 

\mm 

P(15,9) 

j^QQI 

AP(  9,5) 

20.2 

2na 

P(ll|,l|) 

6.6 

P(11.1) 

P(10,5) 

5.6 

5rd 

P(  9,7) 

6.2 

P(  5,7) 

■9 

P(11.7) 

4.5 

tith 

POl.5) 

5.11 

6P(  8,6) 

'  4.4 

R 

1.9 

12 

5th 

6P(  8,2) 

5.4 

P(15,5) 

1.0 

X 

2.5 

ith 

♦ 

a.5 
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product  applicable  to  any  obaervation  time  for  which  aurfaoe  data  are  available.  Accord¬ 
ingly.  one  aoreening-regresaion  run  waa  made  on  the  full  aample  of  720  oaaea  (unatratified) 
in  which  SOO-mb  height,  lOOO-SOO-mb  thiokneaa,  and  the  time  ohangea  of  theae  variablea 
(aee  Table  2-3)  were  eliminated  aa  poaaible  prediotora.  Table  S-6  liata  the  aelected 
predictora,  in  the  order  of  their  aeleotion. 

The  oon4)arative  reaulta  of  the  two  teata.  excluding  and  including  iqpper-air  data, 
are  ahown  in  Tablea  3-7  (dependent  aample)  and  3-8  (independent  aample).  Although 
generally  better  reaulta  were  obtained  when  upper-air  data  were  included,  the  differencea 
do  not  am>ear  to  be  exceaaive,  and  the  uae  cf  regreaaim  equationa  that  exclude  tq)per-air 
data  e;>pears  to  be  of  value  in  aituationa  where  en>cr-air  data  are  unavailable. 

3.2  Aaian  Cyclones 

Aaia  (area  II  of  Fig.  2-1)  was  divided  by  latitude  4(rN  and  longitude  14S*E  into  four 
fairly  even  sones  providing  a  separation  of  the  relatively  youxig,  western  cyclones  from 
the  older,  eastern  cyclones. 

3.2.1  Synoptic  Climatology 

Table  3-9  contains  the  means  and  standard  deviations  of  the  northward  and  east¬ 
ward  displacements  and  changes  in  central  pressure  of  1091  Asian  cyclones,  for  12,  24, 
and  36  hr.  Figure  3-8  depicts  the  mean  cyclone  tracks  in  each  of  the  four  sones.  The 
southern  cyclones  tend  to  move  much  faster  than  the  northern  cyclones  (iqiproximately 
27  knots  in  the  south,  17  knots  in  the  north).  The  maximum  deepening  was  also  noted  in 
the  southern  sones,  with  lesser  deepening  over  the  northern  sones.  In  fact,  after  24  hr 
in  the  northeastern  sone,  there  was  a  slight  tendency  for  filling.  This  is  not  surprising 
inasmuch  as  cyclones  in  the  northeastern  smie  are  usually  in  or  nearing  their  occluded 
phase. 

The  associated  mean  nuq>s  of  pressure,  height  and  thickness  are  shown  as  Figs.  3-9 
through  3-14.  Comparison  of  the  sea-level  pressure  maps  (Fig.  3-9)  reveals  a  cyclone's 
stage  of  development;  i.e.,  the  farther  east  the  cyclone,  the  lower  its  central  pressure. 

The  younger  cyclones,  which  comprise  much  otthe  data  samples  of  the  northwestern  and 
southwestern  zones,  have  mean  central  pressures  of  1004  and  1006  mb,  respectively. 
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TABtE  9-9 

CmiWCTERISTICS  OF  MINTER  CYCLONES  OVER  ASIA 
199J-1999  (0EPW06NT  SANI^LE) 


Zona 

For«eatt 

lnt«rval« 

hr 

ObMTwad 
northaard 
ditplacaawntt 
d^.  lat. 

Obaarvad 

aaataard 
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d^.  lat. 
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mb 
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12 

0.l4 

2.34 
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2.34 

-  2.79 
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NH 

2l« 
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5.55 

-  7.50 

3.60 
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8.34 

36 

1.55 

4.54 
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11.20 
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1.66 
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-  2.95 
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-  2.75 

7.31 

B 

24 
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4.12 

-  5.66 

4.67 

-  3.55 

10.80 

■ 

56 

4.72 

5.56 

-  8.51 

6.70 

-  2.85 

15.74 

12 

1.90 

2.09 

-  4.92 

2.42 

-  5.45 

9.49 

SM 

24 

4.05 

5.32 

-  9.78 

4.i4 

-10.99 

9.26 

56 

6.98 

4.75 

-i4.i6 

5.94 

-i4.87 

11.57 

12 

2.45 

2.58 

-  5.40 

5.06 

-  4.82 

6.59 

SE 

24 

4.79 

5.89 

-10.56 

5.12 

-  8.91 

9.77 

56 

7.45 

5.51 

-i4.58 

7.07 

-12.05 

12.59 

12 

1.55 

2.49 

-  3.97 

2.79 

-  3.61 

6.44 

All 

24 

5.15 

4.10 

-  7.69 

4.85 

-  6.4i 

10.20 

56 

4.77 

5.56 

-11.11 

6.79 

-  8.05 

15.42 

*Negative  valuas  raprasant  daapaning. 
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g.  3-9.  Mean  maps  of  sea-’level  pressure  for  winter  cyclones  in  Asia*  1956-19^  (dependoit  stoiple). 
's  are  labeled  in  millibars. 
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san^le).  Isallobars  are  labeled  In  millibars. 


Fig.  3-11.  Mean  maps  of  500-nd>  height  for  winter  cyclones  in  Asia.  1955-1959  (dependent  sanmle).  1 
are  labeled  in  tens  of  feet. 


The  more  mature  oycloneB,  which  dominate  the  northeastern  and  southeastern  zoues, 
have  central  pressures  of  987  and  996  mb,  respectively.  The  four  12-hr  pressure-change 
charts  (Fig.  3-10)  have  similar  isallobaric  characterlstics-a  katallobarlc  center  about 
one  grid  interval  east  and  an  anallobaric  center  about  four  grid  intervals  west.  At  500  mb 
(Fig.  3-11),  the  height  contours  have  larger  amplitude  in  the  two  northern  zones  than  in 
the  southern.  The  same  can  be  said  for  the  lOOO-SOO-mb  thickness  field  (Fig.  3-13). 

The  absolute  values  of  the  thickness  contours  are  about  1000  ft  lower  in  the  northern  zones 
than  the  southern,  which  is  simply  a  reflection  of  the  temperature  variatitm  with  latitude 
in  the  troposphere.  There  is  similarity  among  the  500-mb  height-change  charts  (Fig.  3-12). 
The  center  of  maximum  falls  ia  near  or  slightly  west  of  the  cyclcme  center.  The  rather 
small  magnitudes  of  both  maxima  and  minima  can  be  accounted  for  only  by  the  apparent 
variability  in  location  of  the  rise  and  fall  centers  among  individual  cases,  tending  to  damp 
them  in  the  mean.  The  12-hr  thickness  change  for  the  1000- 500- mb  layer  (Fig.  3-14)  is 
proportional  to  the  changes  in  mean  tenu>erature  for  the  layer.  Furthermore,  it  is  likely 
that  this  tenu^erature  change  provides  a  rough  estimate  of  thermal  advection  for  the  layer. 
In  the  mean-thickness -change  maps,  negative  values  refer  to  cold  advection  and  positive 
values  to  warm.  Of  particular  interest  is  the  contrast  between  the  cyclones  in  the  south¬ 
western  and  northeastern  zones  of  Fig.  3-14.  The  former,  due  to  its  geographical  location, 
may  be  thought  of  as  a  young  cyclone,  and  the  latter  a  mature  («e.  Note  how  the  axis 
between  the  centers  of  cold  and  warm  advection  shifts  from  an  east-west  orientation  in  the 
southwestern  zone  to  a  northeast-southwest  orientation  in  the  northeastern  zone.  This  is 
consistent  with  life-cycle  cyclone  models,  in  which  the  mature  cyclone  has  polar  air  pene¬ 
trating  equatorward  to  its  rear  and  warm  air  being  transported  poleward  in  advance  of  it. 

3.2.2  Results  of  Using  Both  Surface  and  Upper-air  Predictors 

Regression  equations  were  derived  for  the  four  zones  s^arately  and,  as  in  the 
experiment  with  European  cyclones,  the  imstratified  sample  of  the  entire  area  was  also 
screened.  The  predictands  and  possible  predictors  are  the  same  as  were  used  for  Eunq)e 
(Tables  2-2  and  2-3). 

The  results  of  the  screening  regression  are  listed  in  Tables  3-10  and  3-11. 

Table  3-10  lists  the  selected  predictors,  in  order  of  their  selection  by  the  screening 
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TABLt  MO 
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*lncluding  upp«r-air  predictors 
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procedure.  Md  ttie  jMrocotag*  of  the  total  variance  explained  by  aaoh.  Table  3-11  aum- 
mariaea  the  initial  and  reaidual  atandard  deviatlona,  along  with  the  percent  reduotiona. 

One  intereating  feature  (tf  Table  3-10  ia  that,  for  a  particular  prediotaiid.  there  ia 
a  tendency  for  the  aoreening-regreaaion  program  to  aeleot  the  aame  or  aimilar  predtetora, 
whether  it  be  for  a  lag  of  12.  24.  or  36  hr.  For  exaii4>le.  for  the  northward  component  of 
diaplaoement  (ft)  in  Table  3- 10(e).  the  firat  predictor  aelected  ia  Z(1S.  5)  for  all  three 
time  perioda.  The  aame  holda  true  for  the  eaatward  component  (fi).  for  which  Z(7. 1)  waa 
aeleoted  llrat,  and  central  preaaure  (6).  for  which  ^(6. 3)  waa  aelected  firat.  In  addition, 
aeveral  of  the  aecond  prediotora  aelected  are  either  the  aame  or  aimilar.  Theae  pre- 
dictora  are  like  thoae  aelected  for  the  northern  zone  of  Europe  (aee  Sec.  3.1.2). 

Table  3-12  containa  the  rma  errora  found  in  the  independent  aanqple.  Note  that,  aa 
in  the  caae  of  European  oyolonea.  the  unatrati&ed  equationa  congmte  well  with  the  atratified 
equationa. 

Category  in  teating  on  cyclonea  in  Eurasia  la  being  done  by  an  Air  Weather  Service 
(U8AF)  Weather  Central  in  Japan.  Preliminary  reaulta  on  approximately  100  caaea  for 
the  1962- 1B63  winter  season  are  available  and  are  shown  in  Table  3-13.  The  regression 
equations  for  total  displacement  had  smaller  rms  errors  than  the  Weather  Central's 
operational  forecasts  for  all  three  time  periods,  although  the  36-hr  rms  errors  were 
nearly  equal.  For  central  pressure,  the  regression  method  is  competitive  with  the 
Central's  prediction  through  24  hr. 

3.2.3  Results  of  Using  Surface  Predictors  Only 

A  screening-regressicm  run  was  made  using  surface  data  only,  in  the  same  nuumer 
as  for  the  European  cyclones.  The  predictors,  in  order  of  their  selection,  are  listed  in 
Table  3-14.  Table  3-15  (dependent  sample)  shows  the  percent  reduction  of  variance  and 
residual  standard  deviation.  Table  3-16  (independent  sanqtle)  compares  the  rms  errors 
involved  in  including  and  excluding  upper-air  data.  As  in  the  Eurqpean  test,  better  results 
were  obtained  when  upper-air  predictors  were  included,  but  the  surface-data  regressicm 
equations  appear  to  be  useful  should  u{q>er-air  data  be  unavailable. 
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3.3  European  Anticyclones 

The  test  on  European  anticyclones  reported  herein  represents  an  initial  feasibility 
test  on  anticyclones.  Although  the  anticyclone  has  received  recent  attention  (llj  as  a 
significant  feature  in  a  moving-coordinate  statistical-prediction  model,  the  form  of  the 
predictand  to  use  is  not  so  obvious  as  it  is  in  the  cyclone  problem.  An  anticyclone 
obviously  cannot  be  adequately  described  simply  by  its  central  pressure  and  location. 

The  specification  of  a  predictand  that  describes  the  sluq>e  of  an  anticyclone  or  the  orienta¬ 
tion  of  its  major  axis  could  be  rewarding.  For  this  particular  study,  the  first  attempt, 
however,  was  a  straightforward  application  of  the  technique  to  predict  only  displacement 
and  change  In  central  pressure. 

The  geographical  area  for  which  cases  were  selected  is  area  I  of  Fig.  2-1,  the  same 
as  for  European  cyclones.  In  this  case,  however,  no  attempt  was  made  at  stratification, 
nor  was  there  an  experiment  to  use  surface  predictors  only,  although,  based  on  results  of 
surface-data  regression  for  cyclones  in  Europe  and  Asia,  such  an  £q)proach  seems  entirely 
within  reason. 

3.3.1  Synoptic  Climatology 

Table  3-17  contains  the  means  and  standard  deviations  of  the  northward  and  east¬ 
ward  displacements  and  changes  in  central  pressure  for  658  European  anticyclones  for 


TABLE  5-17 

CmRACTERISTICS  OF  WINTER  ANTICYCLONES  OVER  EUROPE, 
1955-1959  (DEPENDENT  SAUPLE) 


Forecast 
i  nterva 1 , 
hr 

Northward 
displacement, 
deg.  lat. 

Eastward 
displacement, 
deg.  lat. 

Change  in 
centra  1 
pressure,  mb 

Mean 

Std.  dev. 

Mean 

Std.  dev. 

Mean 

Std.  dev. 

12 

-0.44 

2.05 

-1.90 

2.50 

o.4o 

2.57 

24 

-0.91 

3.29 

-5.87 

4.47 

0.59 

4.06 

36 

-1.35 

4.26 

-5.83 

6.57 

0.57 

5.39 
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12,  24,  and  36  hr.  It  can  be  aaen  that  the  mean  track  is  roughly  eaet-southeaetward  at 
about  10  knots,  with  little  ohange  in  oentral  preesure. 

The  mean  maps  of  pressure,  height,  anfl  thickness  are  shown  in  Figs.  3-15  through 
3-20.  As  ejected,  the  features  are  about  reversed  from  those  on  the  cyclone  ma(M.  The 
mean  oentral  pressure  of  the  antioyclones  (Fig.  3-15)  is  1032  mb.  The  12-hr  pressure 
ohange  map  (Fig.  3-16)  shows  a  northwest-southeast  orientation  of  isallobaric  craters, 
with  very  small  magnitudes.  At  500  mb  (Fig.  3-17),  there  is  a  north-south  ridge  about  380 
km  west  of  the  anticyclone.  The  anticyclone  is  located  in  the  northwestward  flow  ahead 
of  the  ridge.  Heiid^t  and  thickness  rises  tend  to  dominate  the  12-hr  change  maps  (Figs. 
3-18  and  3-20). 

3.3.2  Results  of  Using  Both  Surface  and  Upper-air  Predictors 

The  regression  equations  for  European  anticj'clones  were  derived  from  the  entire 
data  sample  (unstratified  and  containing  both  surface  and  upper-air  predictors).  The 
predictors  are  listed  in  Table  2-3.  The  results  of  the  regression  analysis  are  shown  in 
Tables  3-18  and  3-19. 

In  Table  3-18,  the  first  predictor  selected  by  the  screening  procedure  for  the 
northward  component  is  the  same  for  all  three  time  periods,  Z(13,  5).  Its  relatioaship 
with  the  predictand  is  to  favor  greater  southward  displacements  the  lower  the  height  is. 
Here  again,  there  is  evidence  of  a  kind  of  "steering"  predictor.  For  eastward  displace¬ 
ment,  the  central  pressure  was  picked  up  first  for  the  24-  and  36-hr  forecast  intervals, 
and  the  pressure  one  grid  interval  west  was  selected  first  for  the  12-hr  interval.  The 
relationship  is  that  weaker  anticyclones  tend  to  move  eastward  more  rapidly  than  strong 
ones.  For  central  pressure,  there  was  a  different  'leading  predictor"  for  each  of  the 
three  time  periods. 

On  independent  data  (Table  3-20),  the  only  conq>arison  currently  available  is  with 
the  climatology.  The  regressira  equations  show  improvement  over  climatology  for  all 
nine  predictands. 
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4.0  forecast  examples 


Examples  of  forecasts  prepared  for  five  oases  selected  from  the  independent  data 
aanqple  are  given  here  to  acquaint  the  reader  with  the  performance  of  the  various  equa- 
ti<ms  when  iqq>lied  to  some  rather  Interesting  oases.  Representative  values  of  the  over¬ 
all  accuracy  of  the  equations  may  be  found  in  the  verification  statistios  (Sec.  3.0)  based 
on  the  independent-data  test  of  the  entire  independent  sanq>le. 

4.1  Two  European  Cyclones 

Figure  4- 1(a)  shows  a  porticm  of  the  surface  chart  at  initial  time  00  GCT  30  Jan. 

1960.  The  cyclone  of  interest  is  centered  In  extreme-western  U.S.S.R.,  Just  east  of  the 
Baltic  Sea.  The  central  pressure  at  forecast  time  is  1004  mb.  The  SOO-mb  chart  at  the 
initial  time  was  marked  by  moderate  to  strong  sonal  flow  over  central  and  eastern  Europe. 

The  observed  positions  12,  24,  and  36  hr  subsequent  to  forecast  time  (reading 
generally  from  left  to  rlfdit)  are  shown  as  circles  in  Fig.  4- 1(b).  Three  sets  of  12- ,  24-, 
and  36-hr  forecasts  are  also  shown:  triangles  denote  the  forecasts  obtained  by  the  applica¬ 
tion  of  the  equations  derived  from  the  unstratified  develqpmental  sanq>le;  squares  denote 
the  forecasts  based  on  the  stratified  (noxthem-E<me  in  this  instance)  sample;  and  "X's" 
denote  the  forecasts  based  on  equations  cmtaining  only  surface  predictors.  Observed  and 
forecast  central  pressures  are  plotted  adjacent  to  the  corresponding  observed  or  forecast 
position  symbols. 

The  cyclone  under  consideration  here  moved  due  eastward  for  the  first  12  hr  and 
then  moved  east-southeastward  for  the  next  24  hr  of  the  forecast  period.  The  cyclone 
slowly  deepened;  the  central  pressure  had  decreased  to  998  mb  by  the  end  of  the  forecast 
interval. 

The  eastward,  then  east-southeastward  motion  of  the  cyclone  was  forecast  rather 
well  by  the  unstratified  and  surface-predictor  equations  and  the  errors  in  central  pres¬ 
sure  were  rather  small  for  all  three  sets  of  forecasts.  Table  4-1  is  a  summary  of  the 
verification  of  the  18  forecasts  for  this  case. 

Figure  4-2  illustrates  the  forecast  of  an  intense  mature  cyclone.  At  initial  time 
(00  GCT  18  Dec.  1959),  a  966-mb  low  was  centered  near  the  Sietland  Islands,  due  west 
of  Bergen,  Norway.  This  cyclone  moved  nearly  due  eastward  for  the  first  24  hr  of  the 
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forecast  period  and  then  moved  north-northeaetward  at  about  25  knots  for  the  next  12  hr. 
This  cyclone  filled  28  mb  in  36  hr. 

Both  the  stratifled-sample  and  aurfaoe-prediotor  equaticms  yielded  more  accurate 
position  predictions  for  12  and  24  hr  than  the  unstratlfled-sample  equations;  however, 
the  unstratified-sample  equations  gave  the  better  indication  of  northward  change  in  course 
of  the  cyclone  center  toward  the  end  of  the  forecast  period.  All  three  sets  of  equations 
underforecaat  the  24-  and  36-hr  observed  filling,  but  appreciable  filling  was  predicted  by 
all  diree  types  of  equations. 

4.2  Two  Aslan  Cvolwies 

Figure  4-3(a)  illustrates  the  surface  chart  at  initial  time  (12  GCT  5  Jan.  1960)  for 
one  of  two  Asian  cyclone-forecast  examples.  Figure  4-3  (b)  illustrates  the  observed  and 
forecast  tracks  of  the  cyclone. 

At  initial  time,  the  cyclone  was  centered  Just  north  of  Hokkaido,  and  its  central 
pressure  was  984  mb.  During  the  subsequent  36  hr,  the  center  moved  northeastward  to 
the  vicinity  of  the  southern  tip  of  Kamchatka  and  deepened  an  additional  20  mb.  At  500  mb 
at  Initial  time,  there  was  a  noajor  trough  in  the  westerlies  about  10*  of  longitude  west  of 
the  surface  cyclone.  Strong  southwesterly  winds  prevailed  in  the  mid  and  upper  trqpo- 
sphere  over  the  center.  The  verifications  of  the  forecasts  are  given  in  Table  4-3. 

Figure  4-4(a)  shows  the  surface  chart  for  00  GCT  24  Deo.  1959.  Forecasts  of  the 
36-hr  displacement  and  change  in  central  pressure  were  computed  for  the  1002-mb  low 
centered  about  100  naut  mi  east  of  Hc^aido.  As  in  the  previous  example,  the  500-mb 
level  was  mariced  by  a  major  trough  from  10  to  15*  of  longitude  west  of  the  surface  low, 
with  stroig  southwesterly  winds  over  Japan.  In  this  case,  the  SQO-mb  southwesterlies 
divided  into  two  distinct  maxima,  with  the  associated  difiluent  aone  located  over  the 
surface  cyclone.  During  the  ensuing  36  hr,  the  surface  low  moved  northeastward  with  an 
average  speed  of  about  30  knots  and  deepened  45  mb  I 

Figure  4-4(b)  shows  a  plot  of  the  forecast  and  observed  tracks  of  this  rapidly  inten¬ 
sifying  cyclone.  Table  4-4  is  a  tabulaticm  of  the  various  forecast  errors.  Each  of  the  sets 
of  forecasts  gave  good  indications  of  significant  deepening,  and  the  maximum  displacement 
error  for  all  forecasts  was  216  naut  mi. 
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Fi;j,.  First  European  cyclone-forecast  example.  Initial  conditions  at 

ltd  CiCT  .id  Tan.  Cyclone  is  located  at  55.0*N,  22.i?E;  central  pressure  is 

1(111 1  ml). 


60 


X 


Fig.  4- 1(b).  Forecut  verification.  oOburved;  a  uutratified  prediction;  aetret- 
ified  prediction;  xeurfece  predictors  only. 


TABLE  4-1 

VERIFICATION  SWenflY  OF  CYCLONE  SHOWN  IN  FIO.  4-1 


Forecast 

length, 

hr 

Posit  Ion- '/actor  error,  naut  mi 

Cantral -pressure  error,*  ab 

StPOta 

•q». 

Unstrote 

•qta 

Surface 
data  only 

Strat. 

eqa. 

Unstrat. 

•qte 

Surface 
data  only 

12 

102 

120 

36 

-1 

-1 

t3 

24 

228 

132 

60 

0 

-2 

+3 

222 

78 

126 

•2 

0 

t4 

FCantral-prtisur*  arror  aqualb  pradictad  cantrat  prasaura  minus  obsarvad  cantral 
prataura. 


Fig.  4-2(a).  Second  Eurqiean  cyclone-forecast  exaiiq>le.  Initial  conditions  at 
00  GCT  18  Dec.  1959.  Cyclone  is  located  at  61.5*N,  2.5^;  central  pressure  is 
966  mb. 


■iAt<LE  4-2 

VERIFICATION  SUWARY  OF  CYCLONE  SHOWN  IN  FIG.  4-2 


Forecast 

Position-vector  error,  naut  mi 

Central-pressure  error,*  mb  | 

hr 

Strat. 

eqs* 

Unstrate 

aqs. 

Surface 
data  only 

Strafe 

eqs. 

Unstrate 

eqs. 

Surface  j 
data  Oft  I  \  1 

12 

54 

120 

108 

0 

c 

-1 

24 

ia 

106 

90 

-6 

-4 

-7 

56 

216 

180 

254 

-9 

-12 

-11 

•C#ntra I -pressure  error  equals  predicted  central  pressure  minus  observed  centre! 
pressure. 
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Fig.  4-3(a).  First  Asian  cyclone-forecast  example.  Initial  conditions  at  12  GCT 
5  Jan.  1960.  Cyclone  is  located  at  45.0*N,  145.0*E;  central  pressure  is  984  mb. 


Fig.  4-3(b).  ForeoaBt  verification.  oObserved;  eunstratified  prediction;  oatrat- 
ified  prediction;  xaurfaoe  predictors  only. 


TABLE  4-5 

VERIFICATION  SUMMUHY  OF  CYCLONE  SHOW  IN  FIG.  4-3 


36  «  78  24o  -flO  td  fg 


•Cantral-preMure  error  equals  precleted  central  pressure  elnus  observed  centrel 
pressure. 
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24  losa  r  .  .  .  cy^looe-forecast  exainple.  Initial  condition*  at  00  G 

24  Dec.  1960.  Cyclone  is  located  at  42.5*N,  148.(rE;  central  pressure  Is  1002  mb 
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■1  A _ I 


Fig.  4-5(a).  European  anticyclone  forecast  example.  Initial  conditions  at  00  GCT 
19  Mar.  1960.  Anticyclone  is  located  at  63. TN,  12.5*E;  central  pressure  is  1038  mb. 


J 
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Fig.  4-5(b).  Forecast  verification.  oObserved;  nunstratifled  prediction;  ostrat- 
ified  prediction;  xsurface  predictors  only. 


TABLE 

VERIFICATION  SUIWRY  OF  ANTICYCLONE  SHOWN  IN  FIG. 


Forecast 

length, 

hr 

Posit  ion- vector  error,  nout  mi 

Control-pressure  error,*  mb 

Strafe 

•qse 

Unstrat. 

aqs. 

Surface 
data  only 

Strat. 

aqse 

Unctrat, 

eqse 

Surface 
data  only 

12 

- 

66 

- 

- 

-3 

- 

Zk 

- 

102 

- 

- 

-3 

- 

36 

- 

78 

- 

- 

-6 

- 

I  wCantral-pressuro  error  equals  predicted  central  pressure  minus  observed  central 

pressure. 
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4.8  A  European  Anticyoloiie 

Only  one  Mt  of  eguetione  (unetratified  developmental  sample)  waa  derived  for  the 
displacement  and  change  In  central  pressure  of  European  antloyolooes.  Twelve-,  24-. 
and  36-hr  forecasts  are  Illustrated  In  Fig.  4-5(b)  for  the  1038-mb  high  centered  over 
Scandinavia  at  00  GOT  19  Mar.  1960  and  Illustrated  In  Fig.  4-5(a).  This  high-pressure 
cell  drifted  slowly  on  a  rather  erratic  path  for  the  first  24  hr  of  the  forecast  period  and 
then,  during  the  last  12  hr  of  the  period,  moved  somewhat  more  rapidly  southeastward 
and  waa  cmtered  over  western  U.S.8.R.,  Just  east  of  the  Gulf  of  Riga,  at  12  GCT  20  Mar. 
1960. 

The  forecast  track  gave  a  good  Indication  of  the  acceleration  of  the  center  but  under¬ 
forecast  the  Increase  in  central  pressure.  The  vector  errors  of  the  12-,  24-,  and  36-hr 
position  forecasts  were  66,  102.  and  78  naut.  nd.  respectively,  and  the  corresponding 
errors  in  the  central-pressure  forecasts  were  -3,  -3,  and  -6  mb. 
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S.0  SUMMARY  AND  CONCLUSIONS 

Equatlonf  for  the  prediction  of  12-.  24-,  and  36-br  dieplaoeineiite  end  ohaagee  in 
oMitrel  presaure  have  been  derived  for  Aaian  cyclones  and  European  osrclonea  and  aiM- 
oyolonea. 

The  equationa  remained  atable  when  applied  to  an  independent  data  sample. 

Preliminary  oompaiisona  with  subjective  foreoasta  prepared  by  European  and 
Asian  weather  centers  indicate  that  the  equations  yield  results  at  least  as  good  as  those 
obtained  aidsjectively. 

Equations  aiqilicable  to  the  entire  area  omisidered  in  Europe  and  to  the  entire  area 
considered  in  Asia  produced  results  at  least  as  good  as  equations  derived  for  sididivisions 
of  the  areas. 

Equations  were  derived  from  surface  data  only  for  implication  when  upper-air  data 
are  not  available.  Only  a  small  decrease  in  accuracy  was  noted  when  these  equations 
were  conma^^  with  those  derived  from  both  surface  and  upper-air  data. 

The  cyclone-prediction  equations  have  general  applicability  to  cyclones  in  their 
respective  regions.  A  history  of  the  cyclone  track  is  not  required,  nor  are  the  equations 
restricted  to  any  particular  synoptic  pattern. 

^mllcation  of  the  equations  is  conmletely  objective,  and  they  can  be  adapted  for 
use  by  a  weather  detachment*  or  a  completely  automated  weather  central. 

A  useful  by-product  of  the  data  processing  required  for  the  research  is  a  five- 
year  sample  of  error-checked  hemispheric  gridpoint  values  of  surface  pressure  and 
500-mb  heights  filed  in  a  form  suitable  for  direct  electronic  data-processing  machine 
(EDPM)  analysis. 


*See  App.  C  for  worksheets  illustrating  aimlication  at  the  detachment  level. 


6.0  FUTURE  INVESTIGATIONS  AND  EXTENMONS 


Diagnoatio  studlei  should  be  undertaken  to  establish  (pumtitative  relationships 
between  circulation  and  weather.  These  should  coofirm  or  refute  qualitative  relation¬ 
ships  suggested  in  the  literature  and  isolate  those  pressure-pattern  characteristics  rich 
in  weather-predictive  information. 

The  approach  should  be  extended  to  the  investigaticMi  of  the  predictability  of 
Important  characteristics  isolated  in  the  diagnostic  studies  and  should  be  followed  by 
the  development  of  methods  for  combining  the  ensemble  of  significant-feature  predictions 
to  yield  predictions  of  the  entire  surface-pressure  field.  It  is  worth  noting  that  a  major 
portion  of  the  data  required  for  extending  the  research  to  predict  the  surface-pressure 
field  is  now  error-checked  and  available  on  magnetic  tape. 

Investigations  should  be  conducted  to  inoprove  the  equations  presented  in  this  report 
by  defining  new  possible  predictors  with  a  higher  content  of  predictive  information  and 
developing  methods  for  reducing  the  errors  in  estimating  predictor  values. 
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APPENDIX  A.  DATA  PREPROCESSING 


One  of  the  major  probleme  encountered  whm  undertaking  a  itatiatloal  problem  of 
the  nature  described  in  this  report  ccmcerns  the  data  to  be  used  in  the  investigation. 

When  the  scq|>e  of  the  problem  reaches  the  point  at  which  some  11,000.000  gric^int  values 
are  extracted  manually  from  analyzed  nnaps,  tabulated,  and  placed  on  punched  cards,  there 
are  going  to  be  errors  despite  the  moet  careful  and  systematic  error-checking  procedures. 
The  best  that  can  be  hoped  for  is  a  minimization  of  "gross"  errors.  Although  the  data 
received  for  this  project  had  already  undergone  considerable  checks  for  errors,  it  was  felt 
necessary  to  devise  additional  automatic  error-checking  procedures  and  subject  the  data 
to  them.  Although  only  a  small  percentage  of  additional  errors  were  detected,  anmigh  were 
found  to  make  the  procedures  extremely  worth  while. 

There  were  two  basic  data  types  used  in  the  development  of  multiple  linear-regression 
prediction  equations  for  cyclone  and  anticyclone  displacement  and  intensification  over 
Europe  and  Asia.  Basically,  these  are  historical  series  of  weather  maps  (predictor  data) 
and  summaries  of  cyclone  and  anticyclone  tracks  (predictand  data). 

A.  1  Predictor  Data 

The  basic  predictor  data  are  from  the  gridpoint  data  collection  made  for  Project  433 L 
at  the  National  Weather  Records  Center  (NWRC)  at  Asheville.  N.C.  Pressure-height  and 
temperature  at  several  levels  (including  surface,  700  mb,  and  600  mb)  over  the  Northern 
Hemisphere  were  read  from  weather  charts  analyzed  by  the  National  Weather  Analysis 
Center  (NWAC)  on  a  modified  JNWP  grid  of  1020  points.  Each  map  was  actually  read  twice 
to  eliminate  errors.  These  gridpoint  values  were  tabulated,  placed  on  punched  cards,  and 
finally  on  IBM- 705  magnetic  tapes.  The  final  set  of  data  covered  a  continuous  period  of 
twice-daily  (00-  and  12-GCT)  analyses  from  April  1955  to  March  1960,  a  complete  five-year 
sample.  For  the  three  fields  of  interest  in  significant-feature  studies,  sea-level  pressure. 
700-mb  height,  and  500-mb  height,  this  Involved  more  than  11,000,000  numbers.  For  a 
detailed  description  of  the  433 L  gridpoint  data-extraction  procedures  and  card-t^  formats, 
the  reader  is  referred  to  the  Reference  Manual  for  Climatic  Data  Computer  Tapes  [8]. 
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A.  2  Predlctand  Data 

Cyclone  and  anticyclone  tracke  were  tabulated  from  microfilmed  NWAC  surface 
analyses.  The  location  (latitude  and  longitude)  and  central  pressure  of  every  winter 
(November-March)  cyclone  over  Europe  (area  I  of  Fig.  2-1)  was  plotted  at  12-hr  inter- 

I 

vals  for  the  five  winters  corresponding  to  the  grU^int  data.  This  included  cyclones 
that  were  either  previously  or  subsequently  outside  area  I.  Similar  procedures  were 
used  for  Asian  cyclones  and  European  anticyclones.  From  these  plotted  tracks,  only 
cyclones  or  anticyclones  that  were  identifiable  for  36  hr  were  tabulated  and  put  on 
punched  cards.  For  those  that  qualified,  the  following  information  was  tabulated  for 
each  sync^tic  time  (00  or  12  GCT). 

(a)  Date  and  time. 

(b)  Latitude,  longitude,  and  central  pressure  (initial  time  t^). 

(c)  Latitude,  longitude,  and  central  pressure  (t^  +  12  hr). 

(d)  Latitude,  longitude,  and  central  pressure  (t^  +  24  hr). 

(e)  Latitude,  longitude,  and  central  pressure  (t^  +  36  hr). 

The  number  of  cases  shown  in  Table  2-1  is  equivalent  to  the  number  of  cards  that 
were  punched. 

With  the  IBM-705  magnetic  tapes  and  the  punched  cards  of  cyclone  and  anticyclone 
tracks  available  for  processing,  a  variety  of  computer  programs  (Table  A-1)  was  written 
at  the  United  Aircraft  Corporation  Research  Laboratory.*  These  programs  and  their 
sequential  use  are  shown  in  the  flow  diagram.  Fig.  A-1.  In  detail,  the  programs  were 
used  as  follows. 

The  first  step  was  to  process  each  of  the  Asheville  IBM-705  grid  tapes  separately 
by  the  gridpoint  data-conversion  and  -editing  program.  Since  each  tape  contained  only  two 
months  of  grid  data,  there  were  30  tapes  in  all  to  cover  the  full  60-month  (five-yeai)  period. 
This  program  selects  the  variables  to  be  processed:  sea-level  pressure,  700-mb  height, 

*Cosponsored  by  the  United  States  Air  Force,  under  the  433L  Contract,  and  the 
Federal  Aviation  Agency,  under  Contract  FAA/BRD-363. 
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TABLE  A-l 
COMPUTER  PROGRAMS 
WRITTEN  BY 

UNITED  AIRCRAFT  CORPORATION  RESEARCH  LABORATORY 


Program 

Function 

Gridpoint  data 
conversion 
and  editing 

Selects  gridpoint  data  of 

Interest. 

Checks  for  gross  errors. 

Converts  705  tape  to  7090  tape. 

Grid 

correction 

Inserts  corrections. 

Outputs  merged  tape. 

Gridpoint 

Interpolation 

Converts  from  433L  grid  to  JNWP 
grid. 

Outputs  merged  tape. 

Predictand 

preprocessor 

Converts  locations  and  central 
pressure  to  predictand  format. 

Outputs  on  tape. 

S I gn I f 1 cant- feature 
preprocessor 

Generates  a  tape  for  screening 
regression  by 

(a)  selecting  maps, 

(b)  reading  maps  (transforming 
and  interpolating  grids), 

(c)  deriving  predictors. 

Correction* 


Grid 

correction 


Leffcnd 


Q 

O 


Maitnetic  tape 


Punched  cards 


Subroutine 


Separate  profram 


Printed  results 


Q 


706  Aahevllle 
grid  data 


Qrldpolnt 

1 

and  editing 

•V  Edited 
[90 J  grid 
dau 

Merged 

^Wgrld 

1  data 


Cyclone 

traoka 


Corrections 


Qrlt^lnt 

Interpolation 


Predlctand 

preprocessor 


Printed 

errors 


Reordered 
grid  Preprocessed 

data  f  per  )  predlctand  data 

>"  pr^leijS 


Significant- 

feature 

preprocessor 

problem 

Id. 

Screening 

regression 


per  A  Preprocessed 
problem  data 


Prediction 

-equations 


Fig.  A-1.  Computer  programs  used  to  preprocess  data  for  use  In  screening  re¬ 
gression.  Numbers  a|q>earing  within  tq>e  symbols  refer  to  the  number  of  magnetic 
t^s  Involved. 
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and  500-inb  hei^^t.  Other  variables  aa  the  Uqpe,  auoh  as  helf^  at  other  levels  (e.g., 

300  od>)  and  tenqieratures  at  all  levelsi  were  not  processed.  Having  niade  this  selection, 
the  program  then  examined  every  "map"  for  gross  errors.  Differences  in  values  between 
adjacent  gri^iolnta  were  oong>uted  and  a  preselected  percentage  of  the  largest  ones  were 
printed  out,  almig  with  identification.  Also,  grit^ints  that  had  missing  data 

were  noted  as  well  as  values  outside  certain  specified  limits,  and  printed  out.  This  pro¬ 
gram  also  converted  the  selected  data  and  output  it  onto  magnetic  tips  in  a  format 
acceptable  to  the  kind  of  conaputer  to  be  used  for  subsequent  stages  of  processing  (the 
IBM-7090).  The  ouq;>ut  of  this  program  was  a  set  of  three  7090  tapes  for  every  one  that 
was  used  as  liq>ut,  each  of  which  ccmtained  mly  information  from  one  level  (surface,  70(H:id), 
or  500-mb)  for  a  two-m<mth  period.  Thus,  the  processing  of  the  thirty  IBM-70S  tapes 
resulted  in  ninety  IBM- 7090  tapes. 

The  next  step  was  very  inq>ortant  and  was  a  strictly  manual  operation.  It  consisted 
ot  examining  the  printout  of  missing  and  otherwise  suspicious  grldpoint  data.  Miich  at  the 
data,  although  printed  out,  were  acceptable  due  to  unusually  strong  gradients  or  intense 
systems.  For  exang>le,  all  sea-level  pressures  in  excess  of  1060  mb  were  printed  out. 

Many  of  these  were  valid  because  of  the  Siberian  high,  for  instance.  However,  there  were 
some  pressures  between  1090  and  1099  mb  that  were  actually  supposed  to  have  been 
between  990  and  999  mb,  as  shown  by  examination  of  microfilm  m^ps.  The  practice  then 
was  to  go  from  the  printout  to  NWAC  microfflmed  analyses  to  verify  the  validity  of  sus¬ 
picious  data.  In  this  way,  gross  errors  were  noted,  and  corrected  values  were  tabulated. 
Corrections  to  missing  data  were  also  supplied  from  examinatim  of  microfilmed  m^>s. 

It  should  be  enqihasized  that  the  grid  data  received  from  NWRC  were  by  aixl  large  in  very 
good  condition,  with  the  error  percentage  being  extremely  small.  However,  the  nature  of 
the  study  required  even  additional  scrutiny  of  the  data  to  eliminate  so-called  gross  errors. 

The  grldpoint  corrections,  having  been  tabulated,  were  placed  on  punched  cards  to  be 
used  as  input  to  the  grid-correction  program.  All  corrections  cards  for  one  level  for  five 
years  of  <me  month  were  processed  together  (e.g.,  500-mb  heights  for  five  Januaries).  The 
program  inserted  the  corrections  in  the  proper  places  and  wrote  a  new  t^pe  consisting  of 
the  one  level  for  five  years  of  one  month.  Thus  the  program  not  only  Inserted  the  necessary 
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oorreotions  but  also  reduced  the  number  of  required  tapes  from  90  to  36. 

The  grid  format  used  for  the  Initial  data  extraction  at  NWRC  was  a  modified  JNWP 
grid,  in  which  data  were  read  for  every  other  JNWP  gric^int.  It  was  found  desirable  to 
have  data  at  every  JNWP  grit^int  to  facilitate  sid)sequent  computations.  The  grid-inter¬ 
polation  program  acconqilished  this  by  computing  values  at  the  in-between  gridpoints  by  a 
least-squares  fitting  of  a  quadratic  surface,  using  surrounding  data  at  known  gri^^wints. 
Thus,  the  program  output  complete  maps  on  the  1977-point  JNWP  grid.  Another  feature 
of  the  program  is  that  it  used  three  tapes  as  iiq>ut  (e.g.,  five  Januaries  of  sea-level  pres¬ 
sure,  700-mb  height,  snd  500-mb  height).  It  read  in  the  first  record  (map)  from  each  tape, 
which  were  all  for  the  same  time  (e.g.,  00  GCT  1  Jan.  1956),  and  output  them  as  the  first 
three  records  of  the  new  tipe.  Hmice,  the  final  tape  for  a  month  was  in  chronological  order 
(i.e.,  three  levels  for  00  GCT  1  Jan.  1956,  three  levels  for  12  GCT  1  Jan.  1956,  three  levels 
for  00  GCT  2  Jan.  1956,  etc.,  through  31  Jan.  1960).  The  result  was  that  data  on  the  36  tipes 
using  a  modified  JNWP  grid  were  now  contained  on  12  tapes  with  a  full  JNWP  grid,  each 
Upe  representing  a  complete  five-year  month.  These  tipes  were  then  ready  to  serve  as  the 
basic  set  of  "maps"  for  prbdictor  information. 

The  predictand  data  were  handled  independently  of  these  programs.  The  predictand- 
preprocessor  program  accepted  the  punched  card  input  of  significant-feature  tracks  and 
computed  the  two  components  of  displacement  and  the  change  in  central  pressure,  which 
were  the  predictand  forms  desired.  This  informaticm  was  then  placed  on  magnetic  tape. 

An  inoportant  by-product  of  this  program  is  a  printout  of  displacements  and  changes  in 
central  pressure  that  deviate  a  gpreat  deal  from  climatology.  Thus,  an  error-checking 
procedure  was  provided,  and  new  punched  cards  were  prepared  when  errors  were  noted. 

The  predictand  tiq)e  and  the  predictor  tapes  were  then  used  as  iiq>ut  to  the  significant- 
features  yeprocessor  program.  This  program  generated  a  tape  for  screening  regression 
by  selecting  cyclones  (or  anticyclones),  "reading"  m^>s,  and  deriving  predictors.  It  did  so 
in  the  following  manner.  The  first  cyclone  case  was  determined  from  the  cyclone-tabulation 
iiq)ut.  Using  the  time  and  date  of  this  case,  the  gri(^oint-data  tapes  were  searched  to  find 
the  three  data  records  corresponding  to  this  time  (one  record  for  each  vcatp  of  sea-level 
pressure,  700-mb  height,  and  500-mb  height)  and  also  for  the  same  three  maps  12  hr 
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previous  to  that  time.  Theae  alx  records  were  read  into  the  maohlne.  The  next  atq;>  was 
analogous  to  superimposing  a  moving-'coordlnate  grid  (centered  witti  respect  to  the  cyclone) 
over  the  fixed-grid  system  and  interpolating  data  values  for  gritipoints  in  the  movlng- 
ooOTdlnate  system.  The  subroutine  that  computed  these  locations  is  the  grid-transformation 
routine.  Whereas  the  JNWP  or  "fixed"  grid  has  1977  gridpoints,  the  moving  grid  has  only 
221  points.  With  the  use  of  data  at  the  four  surrounding  grid^lnts  of  the  fixed  grid,  values 
were  computed  by  Interpolation  at  each  of  the  221  grit^ints  of  the  moving-grid  system  for 
all  six  maps  stored  in  the  maohlne.  The  procedure  of  reading  maps  into  the  machine,  com¬ 
puting  locations  of  gri4>oints,  and  interpolating  data,  was  acconplished  for  all  cyol<»es  to 
be  considered.  Selection  and  derivation  of  possible  predictors  from  this  newly  generated 
set  of  data  was  then  performed.  Possible  predictors  consisted  of  point  values  of  sea-level 
pressure,  500-mb  height,  1000-500-ntb  thickness,  and  12-hour  changes  of  these  quamtities. 
The  changes  in  latitude  and  longitude  (both  in  degrees  of  latitude)  and  central  pressure  of 
the  cyclone  were  used  as  predictands  for  12,  24,  and  36  hr. 

The  results  of  this  series  of  programs  were  magnetic  t^pes  that  could  be  used 
directly  as  input  to  the  screening-regression  program  to  derive  prediction  equations  for 
cyclones  in  Europe  and  Asia  and  anticyclones  in  Europe.  In  additioa,  useful  by-products 
included  a  relatively  error-free  five-year  set  of  grid^int  data,  which  can  be  used  for 
work  on  significant  features  in  other  seasons,  such  as  cyclones  in  summer,  or  for  other 
meteorological  problems. 


APPENDIX  B.  PREDICTION  EQUATIONS 

The  predlotlon  e<|iiaitl(m«  derived  from  the  regreesion  enalyele  have  the  form 

f  .  A,  4  4  A^X^  4  ...  4  A^X„.  (B-1) 

where  f  ie  the  prediotand,  the  A'e  are  oonatant  ooetfioients  derived  from  the  developmental 
■ample,  and  the  X'a  are  &e  predictors  selected  the  screening  procedure. 

Each  set  of  prediction  equations  consists  of  nine  equations:  the  three  predlctands  of 
northward  displacement,  eastward  displacement,  and  change  in  central  pressure  for  the 
forecast  intervals  of  12,  24.  and  36  hr. 

The  pair  of  numbers  that  is  associated  with  a  given  predictor  in  the  equations  refers 
to  the  grid  location  in  the  (k,i)-grid  system  of  Fig.  2-2.  For  convenience,  the  Sjrmbols  and 
units  of  the  predlctands  and  predictors  are  repeated  in  Table  B-1. 

For  operational  application,  the  equations  in  worksheet  form  represent  a  more  con¬ 
venient  means  for  solution.  An  exanq>le  of  such  worksheets  is  shown  in  Appendix  C. 

B.1  European  Cyclones 

B.1.1  Ectuati<»a  Using  Both  Surface  and  Upper-air  Predictors 
B.1. 1.1  Northern  Zone 

Ni2  -  52.9918  +  0.0129  Z(15.  5)  -  0.0370  Z(9,  7)  -  0.0691P(5, 1)  +  0.0272  Z(13.  S) 

£^2  *  39.5583  -  0.0233  Z(9,l)  +  0.0081Z(13. 9)  -  0.0838  AP(9,  5)  -  0.0178  Z(1S,  1) 

+  0.0268  Z(ll,  7)  -  0.0159  Z(9,  5) 

°12  "  +  0.4022AP(9,5)  -  0.4952P(10,  5)  +  0.3473P(11.  5)  -  0.1160 P(13, 1) 

-  0.2024  AP(7.1)  +  0.0616  AZ(9,  5)  -  0.0800A  0.0377  Z(9, 7) 

^24  *  43.9238  -i-  0.0321Z(16,  5)  -  0.0608Z(9.7)  -  0.0314Z(5, 1)  +  0.0422Z(11, 1) 

-  0.2643^  +  0.0475P(15,7)  -  0.0306H(9,1) 

£34  »  71.2974  -  0.0397  Z(9,l)  +  0.0267  Z(13, 9)  -  0.0340  Z(13, 1)  -  0.0330  AZ(9,  5) 

+  0.0266Z(13.7)  -  0.0515P(10,5)  +  0.0104 Z(5. 11) 

634  “  296.830  +  0.2471  AZ(9, 5)  -  0.4131  P(10, 5)  -  0.1385  AH(9,  5)  -  0.1349X 

•f  0.0652  Z(9, 7) 

£36  “  “  9-7252  +  0.0530Z(15,5)  -  0.0763Z(9,7)  +  0.0270 AZ(13, 3)  -  0.0364Z(5, 1) 

-  0.3310(^  +  0.1294  P(13,l) 
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'  TABLE  B-1 

SYMBOLS  AND  UNITS  USED  IN  PREDICTION  EQUATIONS 
(a)  Pradietanda 


Symbol 

Description 

Unit  of 
meesuremant 

At 

N 

Predicted  northward  displacement 

Deg.  lat. 

t 

Predicted  eastward  displacement 

Oeg.  lat.e 

0 

Predicted  change  In  central  pressure 

mb 

(b)  Pradictors 


Symbol 

Description 

Unit  of 
measurement 

P 

Sea- level  pressure 

mb 

AP 

12-hr  pressure  change 

mb 

Z 

900-mb  height 

10  ft 

AZ 

12-hr  height  change 

10  ft 

H 

1000-;i00-mb  thickness 

10  ft 

12-hr  thickness  change 

10  ft 

Latitude  of  significant  feature 

Deg.  lat. 

A 

Longitude  of  significant  feature 

Oeg.  long.f 

*Negativa  values  Indicate  eastward  displacement. 
tPosItIve  for  W  longitude,  negative  for  E  longitude. 
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-  55.9758  -  0.05912(9.1)  0.0436  2(13.9)  -  0.0396  2(13.1)  -  0.0491  A2(9.  J) 

3o 

+  0.02322(15.7) 

0..  -  252.359  -  0.02892(9.3)  4  0.8375  AP(9. 5)  -  0.7384P(10.  5)  +  0.0893  H(9.  7) 

3o 

-  0.0748  2(13.1)  -I-  0.5178P(11.5) 

B.  1.1.2  Southern  2one 

f),.  -  32.4831  -  0.0637P(7.7)  0.0374  2(13.3)  -  0.0207  2(9.7) 

12 

.  -  74.4198  +  0.0169  2(11.9)  +  0.0755  P(5.  7)  -  0.0187  2(7.1) 

fi,.  -  139.129  +  0.3750AP(9.3)  -  0.3095 P(10.  5)  +  0.1683P(11. 7) 

12 

-  75.8769  -  0.1229 P(7. 5)  +  0.0537  2(13.3)  -  0.0280  2(9.9)  -  0.0350  A2(7. 9) 

24 

*  -  157.067  +  0.03832(11.9)  +  0.1263P(3.7)  -  0.02502(9.5) 

24 

6..  «  286.717  -  0.6027P(10.5)  4  0.4455 AP(9. 3)  *  0.2411P(11. 7)  -  0.1061X 

24 

+  0.0414  H(9, 9) 

N..  ••  132.281  -  0.0439  2(9.9)  *  0.06562(13,3)  >  0.1713P(5.3) 

36 

»  -  323.286  +  0.0490  2(13.9)  0.2245  P(5.  7) 

36 

•  568.932  -  0.6602 P(10.  5)  0.05212(11.9)  -  0.1046X  t  0.4314  AP(9. 3) 

36 

-  0.4810  AP(3. 3) 

B.  1.1.3  Both  Zones 

N..  -  80.4874  -  0.1099P(7.1)  0.0444 P(13, 3)  -  0.1008P(9.7)  0.01982(15.5) 

12 

-  0.0161  H(7, 7)  +  0.0796  P(ll,  3) 

-  53.3866  +  0.02S1P(15.9)  -  0.0638P(11. 1)  +  0.0159P(5.7)  -  0.0673 AP(7.  5) 

12 

-  0.01422(9.1)  +  0.02372(11.7)  -  0.1127P(10.  5)  +  0.1784P(11.  5)  -  0.0165H(13. 3) 
•*■10.00882(5,9)  -  0.1007P(9,  5) 

0^2  >  208.903  +  0.3983  AP(9,  5)  -  0.5119  P(10,  5)  *  0.03702(9.7)  -  0.0894  X 

-*■  0.3955P(11.5)  -  0.1593P(13.3)  ■*■  0.0371  A2(9. 3) 

■  105.753  -  0.1131P(7, 1)  ■*•  0.0294P(15, 3)  -  0.0067P(9,9)  0.02632(15,5) 

-  0.0359  2(9,7)  +  0.04122(11,3)  -  0.0355  A2(9, 9)  -  0.0227H(7.3)  -  0.0898P(11,  7) 
0.0576  P(15,  7) 

-  67.6409  +  0.0313  P(15, 9)  -  0.0763  P(ll.  1)  •)■  0.0492  P(5. 7)  -  0.0253  2(9.1) 

-  0.0062  2(13,9)  -  0.0279  A2(7,  5)  +  0.0240  2(11,9)  -  0.0354  H(13. 3)  +  0.03662(13,7) 

-  0.0630  P(  10,  5) 
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6^^  -  391.755  0.3859  AP(9.  5)  -  0.7520  P(10.  5)  -t-  0.0551  Z(9. 7)  -  0.1249  X 

+  0. 1883 AP(1 1.3)  *  0.4173P(11.5)  -  0.1554P(13. 1)  0.0615AZ(9. 5) 

•t-  0.1603AP(11.9) 

Ngg  -  53.0687  -  0.1641P(7.1)  -  0.0500  P(9, 9)  +  0.0582  Z(15, 6)  -  0.0479  Z(9,  7) 

<«■  0.1319P(13. 1)  -  0.0463  AZ(9. 9)  -I-  0.0420  AZ (13. 3) 

-  -  3.9246  +  0.0414P(15,9)  -  0.0593 P(ll.  1)  +  0.0800 P(5. 7)  -  0.0275Z(9. 1) 

•f  0.0477Z(13.9)  -  0.1078AP(9,5)  -  0.0312Z(13, 1)  -  0.1411P(9.3)  +  0.0993 P(13,  5) 
^36  "  496*676  -  0.6125P(10.  6)  0.5398 AP(9,  5)  *••  0.0227Z(13.9)  >  0.1535X 

+  0.0878  Z(9. 7)  -t-  0.0782  AZ(1 1.3)  -  0.0442  H(ll,  3)  0.0759  AZ(11. 7) 

B.1.2  Equations  Using  Surface  Predictors  Only  (Both  Zones) 

“  37.880  -  0.1313P(7,1)  +  0.0349 P(14, 4)  -  0.0462P(9,7)  +  0.1745P(11.  3) 

4  0.0938AP(8.  2)  -  0.0710^  -  0.1213P(10. 6)  +  0.0206 P(5. 13)  +  0.0335P(17,  7) 

■  -  65.852  +  0.0332P(15,9)  -  0.0872P(il,  1)  +  0.0512P(5,7)  -  0.1282 AP(8, 6) 

+  0.0371P(13,  5)  -  0.1250P(10,5)  +  0.1240 AP(11,  6)  0.1123P(12, 6) 

4  0.0396  P(5,l) 

■  150*26  4  0.3929  AP(9, 5)  -  0.5351P(10,  5)  4  0.1693P(11,  7)  -  0.1612^ 

-  0.0628  X  4  0.2688P(10,4)  4  0.1867  AP(11,  5)  -  0.0498P(3,9) 

^24  "  68.515  -  0.2075P{7, 1)  4  0.0392 P(15,  5)  -  0.0069 P(ll,  9)  -  0.1381^ 

4  0.0561  P(7, 13)  4  0.1251  AP(9,1)  4  0.0989  P(13, 1)  -  0.2152  Pll,  7)  4  0.0914  P(15,  7) 
4  0.1276P(11,3)  -  0.0446P(7,9)  4  0.1138AP(13,  7) 

^24  *  ■  96.546  4  0.0988P(15, 9)  -  0.2067 P(ll,  1)  4  0.0345P(5,  7)  -  0.1538 AP(8, 6) 

-  0.1887  P(10,  5)  4  0.1076 P(ll,  5)  4  0.1566  AP(12, 4)  4  0.1217P(13,  5) 

4  0.0710  P(5.1) 

-  224.82  4  0.5714  AP(9,5)  -  0.6890  P(10,  5)  4  0.0846  P(12,  8)  -  0.1958^»  — 

4  0.3341  P(ll,  5)  -  0.0702  X  4  0.2378  AP(11, 3) 

Sgg  «  93.061  -  0.2107P(7,1)  -  0.0677P(9,9)  4  0.2018P(13.  3)  4  0.0625 P(7, 13) 

-  0.1684<#i  -  0.1962P(11,7)  4  0.1266P(17,  7)  4  0.1936AP(13,  7)  -  0.1426  AP(9, 9) 
figg  -  -  63.968  4  0.1749 P(15, 9)  -  0.2642P(11, 1)  4  0.1113P(5,7)  -  0.3348P(10, 5) 

4  0.3656  P(ll,  5) 
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621.36  -  0.6580  P(10, 5)  +  0.6124  AP(9, 5)  +  0.0090  P(13. 9)  -  0.20060 
+  0.3246P(12,6)  -  0.2081  P(15, 1)  -  0.0904X  -  0.0793 P(5, 13)  +  0.2486 AP(11, 3) 


B.2  Aslan  Cyclones 

B.2.1  Equations  Using  Both  Surface  and  Upper-air  Predictors 
B.2. 1.1  Northwestern  Zone 
=  -  31.4556  +  0.0365  Z(13.  5)  -  0.0192  2(9,7) 

=  -  1.6969  4  0.0052  2(13,9)  -  0.0265H(7, 1)  4  0.02122(11,7) 

=  -  52.9586  4  0.3259  A P(9. 3)  4  0.30650  4  0.0299  2(7,7)  4  0.2310  AP(9,  9)  - 

-  0.0688 P(3,3)  4  0.0748P(15,9)  4  0.1712 AP(13,  5)  4  0.0364  AH(9,  7)  -  0.1964 AP(7, 3 

-  0.3837  P(10,  5)  4  0.3616  P(ll.  5) 

^24  "  ■  *  0.0379  2(15,  .3)  -  0.0224  2(7,9)  4  0.0292  H(13,  7)  4  0.1503  P(ll,  1) 

-  0.1169  P(9,  7)  4  0.0547  P(3. 9) 

=  -  10.7959  4  0.0142  2(13,9)  -  0.0574  H(7,l)  4  0.0489  H(ll,  7)  4  0.0393  AZ(5,  9) 

4  0.1450  AP(11, 7)  -  0.0368  AZ(7,  5) 

=  90.2693  4  0.1809 P(15.9)  -  0.1508P(3,7)  4  0.5431  AP(9, 3)  -  0.1095 H(ll,  1) 

4  0.0593  H(9,  5)  -  0.8853  P(10, 5)  4  0.7945  P(ll,  5)  4  0.0344  2(5,9) 

Ngg  =  -  97.9147  4  0.0522  2(15,3)  -  0.03512(7,9)  4  0.0382  H(13,  7)  4  0.0456  AH(11, 3) 

-  0.1611  AP(9,  9)  -  0.0459  AZ(5,  3) 

Egg  =  -  39.4679  4  0.0305  2(15,9)  -  0.0505  2(7,3)  4  0.03912(11,7) 

Dgg  =  442.266  4  0.2597  P(15.  9)  -  0.2037P(3,7)  -  0.9115P(10,  5)  4  0.5120P(9,3) 

4  0.4123  AP(11,  7)  4  0.3436X  4  0.3879  AP(  11, 3)  4  0.08052(9,7)  -  0.1064  2(5,1) 

B.2. 1.2  Northeastern  Zone 

^12  ”  11.6227  4  0.0363  2(13,5)  -  0.0166  2(9,7)  -  0.0151  AH(5,  3)  -  0.0858P(9,7) 

4  0.0619  P(ll,  3)  -  0.0134  H(7. 9) 

E^g  =  13.8955  -  0.03152(9,1)  4  0.0235  2(11,7) 

^12  218.247  4  0.1788AP(11,3)  -  0.4341  P(10,  5)  4  0.3792  P(ll,  7)  4  0.3394  AP(  11,  5) 

-  0.1664  P(3,  3) 

^24  ^  0.0417  2(13,5)  -  0.0498  2(9,7)  4  0.0423  H(ll,  7)  4  0.1084  P(13, 3) 

-  0.0174  2(9,11)  -  0.0209  H(5,  5)  -  0.14950 
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^24  *  ~  0.0363  2(9,1)  4-  0.06612(13.7)  -  0.060611(13.6)  4  0.1200 P(3. 3) 

-  0.0912P(9.6) 

6^4  "  43.1036  4  0.4399  AP(0. 3)  -  0.0630 H(ll.  1)  4  0.0409  2(11.9)  -  0.7968P(10.  5) 

4  0.4279  P(ll,  3)  -  0.0886  AH(5.  7)  4  0.3330  AP(11.  6)  4  0.2967  P(ll,  7) 

f}..  -  127.871  4  0.0399  2(13.6)  -  0.0061  2(9.7)  4  0.0307  2(15.7)  -  0.1991  P(ll,  9) 

SO 

-  0.0179 H(5.  6)  -  0.0559A2(9.9)  4  0.1168P(13. 5)  -  0.0388AH(9, 5)  -  0.0297H(7,9) 

-  0.03532(3, 1)  -  0.1999^ 

-  -  74.3317  -  0.1068  2(9.1)  4  0.0267  2(13.9)  4  0.1603P(15.  7)  -  0.03122(15,3) 

4  0.0612  H(9, 1) 

-  256.719  4  0.4854  AP(9. 3)  -  0.7068P(10,  5)  4  0.3633  AP(11, 3)  4  0.4300 P(13,  5) 

SO 

-  0.1017  H(ll,  1)  -  0.0104  2(11,9)  4  0.0453  2(15,11)  4  0.0756  2(9.9) 

B.2.1.3  Southwestern  Zone 

«  -  17.3263  -  0.2242  AP(11. 7)  4  0.0315H(13.  7)  -  0.02082(9.7) 

-  53.4754  4  0.01182(11,11)  -  0.0218H(5,7)  4  0.2077  AP(13,  5)  4  0.0347  2(11,7) 

-  0.0535  2(13.1) 

5^2  -  -  67.4397  4  0.1038A2(9.5)  4  0.03802(11.11)  4  0.3936 AP(11, 7) 

«  158.162  -  0.3092  AP(11,  7)  -  0.1891P(11,  7)  4  0.0467  2(15.7)  -  0.0279  H(7.  7) 

-  126.452  4  0.0173  2(15.11)  4  0.0603 P(9, 11)  -  0.0730  2(3.1)  4  0.0735  A2 (11.  7) 

4  0.0307  2(11,9)  -  0.0738  2(15,1) 

-  321.99  8  4  0.1839  AZ(9.  5)  4  0.05312(11,11)  -  0.4428P(15,  7)  -  0.1467  AZ(5,  7) 

4  0.9770  AP(1 1,7)  -  0.4673  P(10,  5)  4  0.4943  P(13, 7) 

Ngg  -  17.4652  -  0.3254P(11,7)  4  0.0621P(15, 7)  -  0.3407 AP(9, 5)  4  0.04682(15,7) 

4  0.1645P(3,7)  -  0.0867  AZ(9,  7) 

E^g  >  265.246  -  0.1202  2(15,1)  4  0.05732(13.11)  -  0.0779  2(3.3) 

>  359.260  4  0.0954  2(11,11)  -  0.5253  P(15,  7)  -  0.3533  AZ(5,  7)  4  0.2859  AZ(7.  7) 

B.  2  ■  1.4  Southeastern  Zone 

■  -  105.598  4  0.1344  P(15, 7)  -  0.1610 P(ll,  7)  4  0.1325 P(13. 3) 

-  35.9840  4  0.0221  P(13, 11)  -  0.0573  2(5,1)  4  0.02812(11,9) 

"  -  3.5849  4  0.5958 AP(9, 3) 
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■  -  302.042  +  0.3830P(16.7)  -  0.4204P(13, 7)  -  0.2458 AP(9. 7)  +  0.2044 AP(13, 7) 

-  0.0409AZ(n.ll)  ♦  0.1164 P(3, 7)  +  0.02262(15,9)  +  0.1816 P(16. 3) 

■  44.9139  +  0.0480  P(  13, 11)  -  0.1035  2(7,1)  +  0.0534  2(11,9) 

«  357.931  +  0.8334  AP(9.3)  -  0.3643P(10,  5)  +  1.1836AP(11, 1) 

Ngg  ■  -  94.7719  +  0.0609  2(15,7)  -  0.0309  2(11,9)  -  0.2741  P(ll,  7)  +  0.3170 P<15, 3) 

Egg  «  65.7869  +  0.0394 P(  13. 11)  -  0.17572(7,1)  +  0.07712(11,9)  +  0.3502P(13,  5) 

-  0.2724  P(9, 7) 

°36  “  551.109  +  0.7557  AP(9. 3)  -  0.6709  P(10,  5)  +  1.4226  AP(11, 1)  +  0.0623  2(13,9) 

B.2.1.5  All  Zones 

Ni2  =  -  22.7440  +  0.00312(15,5)  -  0.00512(7,7)  +  0.02752(13,5)  -  0.1019 P(ll, 7) 

+  0.0362P(13,  5)  +  0.0375 AP(13, 1)  -  0.0528 AP(9, 7)  +  0.0257P(15,9)  -  0.01362(9,7) 
+  0.1010P(11,3)  -  0.0413  P(10, 5)  -  0.0092  2(5,3)  -  0.0541«^ 

E^2  =  -  48.4427  -  0.0118  2(7,1)  +  0.0139  2(11,9)  +  0.0874  P(13.  5)  -  0.0147  2(11,1) 

-  0.1408P(9,  5)  +  0.0970  AP(  11,  5)  +  0.0536  P(5, 3)  -  0.0160  AZ(7,  7)  +  0.0578 P(ll,  7) 

+  0.0059  2(15, 11) 

Dj2  =  125.717  +  0.1680AP(9,3)  +  0.2612  AP(11,  3)  -  0.5793P(10,  5)  +  0.4346 P(ll,  5) 

+  0.1997AP(9,  5)  -  0.2778  AP(5, 3)  +  0.0955  X  +  0.0599  AZ(9,  7)  +  0.0169  2(11,11) 

+  0.1281AP(13,  5)  -  0.0423  AH(5,  5)  -  0.0371  A2(13,  7) 

^24  ""  13.1651  +  0.02152(15,5)  -  0.0170  2(7,7)  +  0.0218  2(11,3)  -  0.1677  P(ll,  7) 

+  0.1289 P(13,  5)  -  0.01862(5,3)  -  0.0312  A2(9,  7)  +  0.0186  2(15,7)  -  0.0116  2(9,11) 
^24  *  "  ■  0.01912(7,1)  -  0.0049  2(13,9)  +  0.0958  P(13,  7)  -  0.0183  2(9,3) 

+  0.03242(11,9)  +  0.0143  2(15,11)  -  0.0228 A2(7,  5)  -  0.01932(13,1)  -  0.0660P(9,7) 
+  0.0882  P(5,  5)  +  0.1150AP(11,5)  -  0.1289  P(9,  5)  +  0.0935  P(13,  5) 

°24  "  651.448  +  0.5199  AP(9, 3)  -  0.1270H(11,  1)  +  0.0637  H(9,  7)  -  0.8334 P(10,  5) 

+  0.1907  P(ll,  5)  +  0.2060  AP(9,  7)  +  0.1144  \  +  0.1316  AP(11,  3)  +  0.03412(15,11) 

-  0.0542  AH(5,  7)  -  0.3040  AP(5, 3)  +  0.0828  2(11,3)  -  0.1502  H(ll,  7)  +  0.09982(11,7) 
4  0,2869  AP(11,  5) 

N36  =  49.6274  +  0.04712(15,5)  -  0.0205  2(7,9)  -  0.1116  P(ll,  9)  -  0.0124 H(5, 5) 

+  0.1872P(13,5)  -  0.0477  A2(9.  7)  -  0.1220  P(ll,  7)  +  0.0264  2(15,9)  -  0.1506^ 

-  0.0207  2(5,3)  -  0.0167  2(15,11) 


89 


-  85.2273  -  0.0239  2(7.1)  -  0.00282(13.9)  -  0.0472H(7,S)  *  0.1146P(15.7) 

So 

-  0.04812(13.1)  >  0.1860P(9.5)  +  0.04402(11.7)  -t-  0.03102(13,11)  -I-  0.01942(3,7) 

■  470.246  -I-  0.7513^(9.3)  -  0.7642P(10,  5)  -  0.0975 H(ll.  1)  +  0.08322(9.7) 

*  0.0587  2(15.11)  0.4251P(11. 3)  0.2082  X  -  0.2011P(5.5)  +  0.3159^ 

-  0.0618 AH(5. 7) 

B.2.2  Eouatlong  Using  Surface  Predictors  Only  (All  2ones) 

Ni2  -  -  97.697  +  0.0686  P(1S.  S)  -  0.0416  P(6, 7)  -  0.0848<)  4  0.0283  P(17, 9)  -  0.0276 X 

4  0.1276P(12,4)  4  0.0633P(1.3)  -  0.2061P(10. 6)  4  0.1119 P(9,  5)  4  0.0667  AP(10, 2 

-  0.0551  P(7,l) 

fejjj  -  >  16.153  4  0.0513 P(13, 11)  4  0.0818^  -  0.1183P(9, 1)  4  0.0738 AP(11, 3) 

-  0.0579  AP(8,  6)  4  0.1327  P(12, 6)  -  0.1443  P(9, 5)  4  0.0870  P(6, 2)  4  0.0904  AP(11, 1 

-  78.857  4  0.3324AP(9, 3)  4  0.2617  AP(12. 4)  -  0.7230P(10,  5)  4  0.2674P(10, 6) 

-  0.3066  AP(6,  2)  4  0.0661X  4  0.3808P(11,  5)  4  0.1942  AP(10,  5)  4  0.1230  AP(3, 9) 
N24  ■  -  171.36  4  0.1917P(15,5)  -  0.1505P(10,  8)  4  0.0623P(17, 9)  -  0.1088  AP(9,  7) 

-  0.1145<<»  4  0.0657  P(l,  7)  -  0.0477  X  4  0.0913  AP(13, 3) 

^24  “  "  34.132  4  0.0769 P(13, 11)  4  0.1788^  -  0.1993P(11. 1)  4  0.1447P(13. 7) 

4  0.1444AP(12,2)  4  0.0430 P(17, 11)  -  0.1527P(9,5)  4  0.1373 AP(12,  6)  4  0.1053P| 

-  0.1184^(7,5) 

-  275.08  4  0.5172 AP(9, 3)  4  0.1573«^  -  0.1775P(3,3)  4  0.08,34 P(15, 11)  -  0.8181F 
4  0.3756P(10,6)  4  0.4929 AP(12, 4)  4  0.1132X  -  0.3627 AP(6. 2)  4  0.1972 AP(3, 9) 

4  0.2643  P(12. 6)  4  0.2361  AP(9,  5) 

-  -  237.67  4  0.0367P(17,7)  -  0.2026 P(10,  8)  4  0.2274P(15,  5)  -  0.2013(#»  4  0.10031 

-  0.0651  X  4  0.0754  P(17. 9)  4  0.1184  AP(14, 4) 

figg  =  78.152  4  0.1242P(15,11)  4  0.2854<^  -  0.2789P(11, 1)  4  0.1929 P(13,  7) 

4  0.1747 AP(12, 2)  -  0.2679 P(9, 3)  4  0.1292P(5,3)  4  0.1463 AP(12, 6) 

®36  “  434.10  4  0.6815AP(9,3)  -  0.9251P(10,  5)  4  0.1705P(15, 11)  4  0.5924  AP(12, 2) 

4  0.4921  P(12, 6)  4  0.3894  AP(10. 6)  -  0.4649  AP(5, 3)  4  0.1711X  -  0.1558  P(5, 9) 

4  0.2665  AP(5, 9) 
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B.8  luronaia  Antlcvelone« 

B«8.1  BauitloM  UaiiMt  Bofli  Burttoe  and  Upper-tlr  Predlotori  (Both  Zones) 

^18  ”  *  0.02182(18.5)  *  0.0740  AP(9, 7)  +  0.0746  AP(11.  7)  -  0.0129  2(9,3) 

*  0.0091  2(11.9)  -  0.18S4AP(9.8)  0.0158A2(9.9)  -  0.0405^(5.5)  -  0.0105H(11.  7) 

*  0.0121X  >  0.0197P(9. 11)  -  0.1114P(11.3)  +  0.1059P(11. 7) 

«  -  85.2348  *  0.0837P(9.5)  0.1579  AP(7. 7)  -  0.0979  AP(13.  7)  -  0.02162(13,1) 

*  0.0211  H(7. 7) 

•  122.741  *  0.2811AP(9.5)  -i-  0.1235 AP(11, 7)  -  0.4506P(10, 5)  0.0695^ 

•t-  0.3481P(11.5)  *  0.0449AH(9.S)  0.015842(15,7)  +  0.12084P(9, 7)  -  0.0097H(13,3) 

-  13.2920  4  0.0251  2(13.5)  4  0.21134P(9. 7)  4  0.10044P(11. 9)  4  0.0434 P(15, 9) 

-  0.0894P(11.1)  -  0.00712(5,11)  -  0.1539 4P(11, 3) 

t.  *  -  70.3194  4  0.1549 P(10,  5)  4  0.2875 4P(7,  5)  -  0.16614P(11.  5)  r  0.1028P(5.5) 

24 

-  0.10804P(13.9)  4  0.01122(11,9)  -  0.0332  2(13,1)  4  0.02892(7,5)  4  0.16194P(7.  7) 

■  201.067  4  0.2287 AP(11, 7)  -  0.01542(9,5)  4  0.22624P(9, 7)  -  0.6209 P(10, 5) 

4  0.4493  P(ll,  5)  4  0.1280^  4  0.3399  4P(9, 5) 

■■  75.5539  4  0.0449  2(13,5)  4  0.2304 AP(9, 7)  -  0.1739 P(ll,  1)  -  0.00352(5,11) 

36 

4  0.1337P(11,9)  -  0.1170AP(5,5)  >  0.0760 P(9, 11)  >  0.0429 A2(9, 1) 

4  0.0351^2(9,9)  -  0.0191  H(7, 7) 

E,-  ■  -  133.905  4  0.1773 P(10,  5)  4  0.3652^(7,5)  4  0.0581X  -  0.1780  AP(11,  7) 

36 

-  0.0338  2(13,1)  4  0.02272(13,11)  -  0.0500P(3,7)  4  0.0461  2(9,7)  -  0.0344  2(9,1) 

4  0.0542  A2(7, 7) 

^36  ■  359.208  -  0.02102(9.7)  4  0.6058  AP(9, 7)  -  0.3184P(10,  5)  4  0.1777^ 

4  0.2401 4P(13. 7) 

B.3.2  Equations  Using  Surface  Predictors  Only  (Both  2<»ea> 

Not  conq>uted. 
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APPENDIX  C.  PROCEmjRE  FOR  OPERATIONAL  PREMCTION 


C.  1  CoMtructlon  of  Wortoheett 

The  following  exani|>le  demonstrates  the  procedure  for  making  an  operational  fore> 
oast  of  the  displacement  and  change  in  central  pressure  of  European  cyclones  for  12,  24, 
and  36  hr.  The  case  to  be  illustrated  is  hypothetical,  but  should  give  a  general  idea  of  how 
worksheets  for  other  equations  can  be  derived. 

For  this  exanq;)le,  the  unstratified  European  set  of  equations  was  chosen  (see  Sec. 

B.  1.1.3).  In  addition  to  arranging  the  predictors  on  the  worksheets  in  a  form  convenient  for 
tabulation  and  conqiutation,  the  following  types  of  modificatitms  are  desirable  in  preparing 
worksheets. 

(a)  Transform  the  predictands  so  that  the  final  result  of  the  conq>utation  yields  the 
forecast  latitude.  Imigitude.  and  central  pressure.  This  can  be  done  easily  by  adding  the 
initial  latitude  and  central  pressure  to  the  N-  and  fi-equations  cf  A^q^ndix  B,  which  in 
effect  is  making  these  initial  conditions  predictors  multiplied  by  a  coefficient  of  1.  If  an 
initial  condition  is  already  a  predictor,  its  coefficient  is  modified  by  adding  1.000  to  it. 

Thus  for  6^2  (12-hr  predictim  of  central  pressure),  P(10, 5),  which  has  a  coefficient  of 
-0.5119  in  Sec.  B.  1.1.3,  becomes  -^>.4881.  The  longitude  is  a  little  more  involved.  There, 
it  is  necessary  to  convert  the  result  from  degrees  of  latitude  to  degrees  of  longitude  before 
adding  on  the  initial  longitude.  Thus  the  final  longitude  can  be  e;q>ressed  by 

Forecast  longitude  ■  X  +  £  sec  - - - ,  (C-1) 

£ 

where  \  is  the  initial  longitude  in  degrees  of  longitude  (negative  if  east  longitude),  £  is  the 

predicted  eastward  displacement  in  degrees  of  latitude,  ^  is  the  initial  latitude  in  degrees 

of  latitude,  and  0  is  the  forecast  latitude  in  degrees  of  latitude, 
r 

(b)  Transform  SOO-mb  height  and  1000-500-mb  thickness  predictors  by  eliminating 
the  ten-thousands  digit  from  them.  This  leaves  a  three-digit  number  (thousands,  hundreds, 
and  tens  of  feet),  which  facilitates  computations.  This  transformation  will  cause  the  con¬ 
stant  additive  term  of  an  equation  involving  height  or  thickness  to  change.  Note  that  no 
such  transformations  are  necessary  for  height  or  thickness  changes. 
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(o)  When  a  1000- 500- mb  thickness  chart  Is  unavailable,  the  predictor  can  be  com¬ 
puted  accurately  enou^  by  converting  the  sea-level  pressure  to  1000-mb  height  by  use 
of  the  hypsometric  formula.  If  a  mean  temperature  of  (fC  between  the  surface  and 
1000  mb  is  assumed  (which  is  reasonable  for  winter  situations),  the  following  equation 
can  be  used: 

H  -  Z  -  2.62(P  -  1000),  (C-2) 

where  H  Is  the  lOOO-SOO-mb  thickness  in  tens  of  feet,  Z  is  the  500-mb  height  in  tens  of 
feet,  and  P  is  the  sea-level  pressure  in  millibars. 

Construction  of  a  grid  overlay  similar  to  that  in  Fig.  2-2  may  be  facilitated  by 
recalling  that  the  grid  has  the  same  spacing  between  points  as  the  JNWP  g^id,  i.e.,  1  in. 
on  a  polar-stereographic  projection,  with  standard  parallel  at  60*N  for  a  map  scale  of 
1:15,000,000. 

C.2  Use  of  Worksheets 

..  C.  2. 1  Data-tabulatton  Worksheet 

Using  the  grid  overlay  to  locate  the  appropriate  points,  enter  the  sea-level  pressures 
and  500-mb  heights  from  the  current  charts  (in  the  correct  units  and  format)  under  t^. 

[To  position  the  grid  overlay,  place  the  point  (k,  i)  *  (10, 5)  over  the  t^-position  of  the 
cyclone,  and  orient  so  that  the  line  k  •>  10  coincides  with  the  meridian  passing  through  the 
tg-position.  ] 

Determine  the  t  entries  for  the  required  points  from  the  12-hr-old  charts. 

[Maintain  the  gridpoint  (k,  I)  >  (10,  5)  over  the  t^-posiUon  of  the  cyclone;  the  values  recorded 
for  t_^2  ^  obtain  time  changes  of  heights  and  pressures  at  the  several  points.  1 

Determine  the  necessary  A's. 

Determine  the  1000-500-mb  thicknesses  at  the  required  points  either  by  using  the 
formulae  given  at  the  bottom  of  the  page  or  by  reading  the  values  directly  from  the  analyzed 
1000-500-mb  thickness  chart  (using  the  grid  as  sm  overlay)  or  from  a  nomogram. 

C.2. 2  Prediction  Worksheets 

Enter  on  each  of  the  three  prediction  worksheets  the  {q>pr(g>riate  parameters  from  the 
data-tabulation  worksheet. 
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Fig.  C-1.  Data-tabulatlon  worksheet. 
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Fig.  C-2.  Twelve-hour  prediction  worksheet. 


96 


WMmM  enrcuMi  rbuotzom  {ak  mat) 


UTXTUDI 

Ut(0) 

ALLi 

1.0000* 

P(  7,1) 

JAikx 

^.1131  • 

K15,3) 

0.029a  - 

P{11,7) 

ioia  X 

.0.0S92  « 

P(15.7) 

JftlLt 

0.0)74  - 

P(  9,9) 

JAllx 

.X).0067  • 

£(11,3) 

M7  X 

0.(V.t2  - 

£(i).S) 

0.0263  • 

£(  9,7) 

JZiLx 

^.0359  « 

U(  9,9) 

3  X 

^.0)))  • 

H(  7,3) 

JlSt.x 

..0.0227  - 

:iUN  OK  «00UCTS  _____ 
C  K3TANT  AOOtnVE  114.6$3 
TCSUCk:^  UTITOUi  .HSLl. 


uwnuni 

P(11,1) 

jam  X 

^.0763  •  _ 

P(10,)) 

JAftLX 

^.0630  -  _ 

P(  ),7) 

Jfiii.X 

0.0a92  -  _ 

P(15,9)  . 

Ufijlix 

0.0313  “  _ 

z(  9,1) 

.JkLx 

-0.02)3  -  - 

2(13,7) 

_tfil  X 

0.0366  -  _ 

2(11,9) 

0.Q2A0  -  _ 

2(13.9) 

7t>  X 

-0.0062  •  _ 

6A(  7,5) 

_ a.1 

-0.0279  -  _ 

H(13,3) 

.Ufi.x 

-0.035A  -  _ 

SUM  cr  psorucTs 


P(»3,1) 

CAKTRAL  PRASS'JjU 
1013  X  .U.1)}A 

P(10,5) 

lOOl  X 

0.2ASC 

P(11,5) 

1007  X 

0.4173 

6P(11,3) 

X 

O.I8S3 

AP{  9,5) 

Q  X 

0.3*59 

4P(11,9) 

Q  X 

0.1603 

Z(  9,7) 

Tl  X 

0.0551 

AZ(  9,5) 

-g  X 

0.0615 

Ion;  (a) 

zJJJL  X 

•0.1249 

SUM  OP  P.iOObCTS 


CCNSTAKT  aOOITTO  A<6.853 


CCMSTANT  AODinVB  6t  .3409 

(1)  Toul __=*L£ia _ 

(2)  1/2  (irltUl  *  f«tt  Ut.)  3^.4 

(3)  1.8  (2)  /.1S7  X  (1)  »  -S.S 

(a)  Inltlkl  Lorcltwl*  (a)  “1 3. 9 

(5)  PCTiECAdT  U)IOITliaB  [(3WA)]’li2 


POST  CENTRAL  FWS'jUW  IOOX  ' 


PUUCAST 

CLiprr^Ticnia  or  sUAixcTiTf 

VSSiPICAnoR 

lAtituO* 

MO.l 

Loi^ltud* 

-iS.1 

Cant.  Pratturt 

IflQA _ 

Fig.  C-3.  Twenty-four-hour  prediction  worksheet. 
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mOKAM  CTOUMC  lUDIOnON  (3i  NOOt) 


Ut(^) 
P(  7,1) 
K13,1) 
K  9,9) 
Z(15,$) 
Z(  9,7) 
AZ(13,3) 
7,9) 


unTuic 

3^.1  X  1 .0000  -  . 

ia/<i  Z  -0.1541  - 
iQlff  Z  0,1319  - 
JQI9  Z  -Q.QMQ* 

Z  0.05M  -  _____ 
7fl  Z  -0.0479  -  _____ 
-3  Z  0.C420  •  _____ 

_ iSLZ  .0.0453- _ 

3UM  Cr  PRODUCTS  _____ 
CONSTANT  AOUITin  63.35R7 
PORKCASr  UTIIUIS  We.v 

GANTRAL  PRBS9IM 


P(10,5)  <OQ»  Z  0.3#75  - _ 

4P(  9,5)  Q  Z  0.5398- _ 

Z(  9,7)  Ttl  X  0.C878  - _ 

2(13,9)  Ttt  Z  0.0227  - _ 

42(11 .3)  X  0.0782  - _ 

42(11,7)  -a.  X  0.0759  - _ 

H(11,3)  _Ifll.Z  .0.0442- _ 

Lone  (4)  -13  *  X  -0.1535  - _ 

SUN  or  PRODUCTS  _____ 
CONSTANT  AODTTIVI  562.976 
POST  CiiMKAL  nUSSURK  1003 


lOQl  X  0.3875  - 

Q  X  0.5398  - 

Ttl  X  0.C878  - 

jui  Z  0.0227  -  . 

-g  X  0.0782  - 

-a.  X  0.0759  - 

X  .0.0442  -  . 


UZnZTUDI 

P(ii,t)  JfliiLX  -0.0593  - _ 

p(  9,3i)  1QI1  X  -0.1411  - _ 

P(13,5)  laia  Z  0.0993  -  — — 

P(  5,7)  lftl%  X  0,0800  - _ 

P(15,9)  loai  X  0.0414  - _ 

4P(  9,5)  Q  X  -0.1078  - _ 

Z{  9,1)  >6*  Z  -0.0275  - _ 

2(13.1)  _1I£.X  -0.031a  -  _____ 

2(13,9)  Z  0.0477  - _ 

SIM  or  PRODUCTS  _____ 
COMSTAMT  AODITm  -14.9246 

(1 )  T«rt*i  -7,;zai _ 

(2)  1/2  (iniUsl  *  feat  Ut.)  37.r 

(3)  nac  (2)  JJLiI_X  (1  )-_=!!_ 

(4)  Initial  Longltiida  M  -13. i 

(5)  PORXCAST  uxxaiw  [(3>t-(4)]':2dka 


PORSaST 


CONVSNTIiiNAL  CR  SUUECTm  I  VBIPICATIOR 


latituda _ 

Lonfltvda 
Cant.  Praawra 


Fig.  C-4.  Thlrty-six-hour  predictiem  worksheet. 
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To  determine  the  foreout  latltudec: 


(a)  Compute  the  products  of  the  Individual  parameters  and  their  respective 
coefficients. 

(b)  Sum  these  products.  (This  can  be  done  simultaneously  with  oonqfniting  the 
products  by  using  the  "Accumulide  Multiply"  key  on  most  desk  calculators.) 

(c)  Add  the  'Vsonstant  additive"  to  this  sum;  the  total  thus  derived  is,  in  each  case, 
the  forecast  latitude. 

To  determine  the  forecast  longitudes: 

(a)  Conq>ute  the  products  of  the  individual  parameters  and  their  respective 
coefficients. 

(b)  Sum  these  products. 

(c)  Add  the  "constant  additive"  to  this  sum;  enter  the  total  thus  derived  on  line  1  on 
the  right-hand  side  of  the  page.  (This  is  the  forecast  westward  {*)  or  eastward  movement 
(-)  in  degrees  of  latitude.) 

(d)  To  convert  the  forecast  movement  to  degrees  of  longitude,  take  one-half  of  the 
sum  of  the  initial  latitude  (^)  and  the  forecast  latitude  (determined  above).  This  gives  a 
mean  latitude. 

(e)  Using  a  secant  table,  determine  the  secant  ("sec")  of  this  mean  latitude.  Enter 
the  value  on  line  3  on  the  right-hand  side  of  the  page  in  the  blank  following  the  words  "sec 
(2)." 

(f)  Multiply  this  secant  value  by  the  value  entered  on  line  1,  i.e.,  by  the  total  cited 
in  (c),  above.  Enter  this  product  at  the  far  right  of  line  3. 

(g)  Enter  on  line  4  the  initial  longitude  of  the  cyclone  at  time  t^.  Note  that  east 
longitude  is  entered  with  a  minus  sign,  west  longitude  with  a  plus. 

(h)  Algebraically  add  together  the  values  on  lines  3  and  4,  and  enter  this  total  on  line 
5;  this  is  the  forecast  longitude.  (Thus,  the  value  on  line  3  is  a  forecast  loi^tude  ohsnge; 
when  added  to  the  original  longitude,  the  forecast  longitude  is  obtained). 
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To  forecast  the  central  pressure  of  the  cyclone,  use  the  same  procedure,  but  with 
values  of  pressure. 

Unusual  displacements  or  changes  in  pressure  should  be  recomputed. 

Enter  the  forecast  latitude,  longitude,  and  central  pressure,  the  conventicmal  or 
subjective  forecasts,  and  the  verification  in  the  box  at  the  bottom  of  the  page. 
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